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We review key experimental and theoretical results on the Barkhausen efltect, focusing on the 
statistical analysis of the noise. We discuss the experimental methods and the material used and 
review recent measurements. The picture emerging from the experimental data is that Barkhausen 
avalanche distributions and power spectra can be described by scaling laws as in critical phenomena. 
In addition, there is growing evidence that soft ferromagnetic bulk materials can be grouped in 
different classes according to the exponent values. Soft thin films still remain to be fully explored 
both experimentally and theoretically. Reviewing theories and models proposed in the recent past 
to account for the scaling properties of the Barkhausen noise, we conclude that the domain wall 
depinning scenario successfully explains most experimental data. Finally, we report a translation 
from German of the original paper by H. Barkhausen. 
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I. INTRODUCTION 



"Eisen gibt beim Unmagnetisionen ein Gerausch; bei ganz gleichmafiigen Anderung der magnetomorischen Kraft 
klappen die Molekularmagnete sprungweise in ihre neue Lage und erzeungen dadurch in einer dariiber geschobenen 
Spule unregelmafiige Induktionsstofie, die sich im Telephon als Gerausch erkenntlich machen" jlj [As the magnetomo- 
tive force is smoothly varied, the molecular magnets flip in jumps to their new position. Because of this, they generate 
irregular induction pulses in a coil wound around the sample, that can then be heard as a noise in a telephone.] (See 
the Appendix for a complete translation of the original paper). 

All the essential elements of the magnetization noise as discovered by H. Barkhausen are contained in his words 
introducing the famous paper of 1919. A piece of iron, the simplest magnetic material, produces a jerky noise when 
magnetized by a field smoothly changing in time as, for instance, given by the slow movement of a magnetic yoke 
driven by a hand. This noise is irregular, in contrast with the regularity of the applied field, and can be easily heard 
with a microphone. This was the first indirect evidence of the existence of magnetic domains, postulated a few years 
before by Weiss 0: actually, Barkhausen believed that the noise was produced by the sudden reversal of the entire 
magnetic domains (the "molecular magnets"). Only several years later, in 1938, Elmore made the first observation 
of the motion of domain boundaries in a cobalt crystal (described theoretically in 1932 by Bloch and thus referred as 
Bloch walls) under the application of an external field. Curiously, Elmore did not recognize this motion as the source 
of the noise heard by Barkhausen. This instead was successively evidenced by Williams and Shockley measuring a 
SiFe single crystal 4|. It is worth noting that in this very year, Kittel in his famous Physical theory of ferromagnetic 
domains, warned "it may be well to correct the widespread notion" that the Barkhausen effect is connected to the 
complete domain reversal, as Williams and Shockley showed "quite clearly [...] that the Barkhausen discontinuities 
correspond to irregular fluctuations in the motion of a domain boundary (Bloch wall)" [from |5l|, pg. 9]. 

Once the origin of Barkhausen (BK) noise was understood, it was soon realized that it could be used as an effective 
probe to investigate and understand the magnetization dynamics in soft magnetic materials, and explain some of 
their hysteresis properties. At the same time, it was also clear that the task was far from being easy. As significantly 
noted "the Barkhausen effect research is both simple and hard" [Ref. @, pg. 46], in the sense that it is relatively 
easy to perform a measurement, but it is much more difficult to interpret it. This is due to the stochastic character 
of domain wall (DW) motion, which proceeds in jumps, or avalanches, when slowly magnetized. This dynamics is 
strongly affected by material microstructure, and other effects such as the demagnetizing field, external stress, etc. 
For this reason, earlier studies were focused on materials having the simplest domain structure, like the SiFe single 
crystals (see for instance 0), where it is possible to obtain a single DW in frame- type samples. 

The intense activity of Barkhausen studies up to the 70s, both from the theoretical and experimental side, are 
extensively described in the reviews of Rudyak Q , and of McClure and Schoeder . Both papers report a large series 
of experimental data on the statistical properties of the noise: distribution of duration and size of the avalanches, 
power spectra, pulse propagation along the sample, etc. The approach followed to interpret the data is fully phe- 
nomenological, without a precise connection with the DW dynamics. For instance, power spectra llOj and cross 
spectra using two separated coils are described by a proper superposition of elementary independent jumps, whose 
spectral characteristic are defined at priori. These elementary jumps are not strictly related to well specified DW 
displacements. Another approach tried to mimic the active interactions on a DW considering the effect of a series of 
springs on a rigid wall, and calculating the resulting power spectrum Cilia. 

Both reviews also present an exhaustive discussion on the possible effects of the measuring setup on the detected 
signal, considering the effect of the dimensions of both the solenoid and the pick-up coils, the different time constants 
related to the electrical circuits, and, more importantly, the effect of the eddy currents. This detailed analysis testifies 
the accuracy and capability in doing BK experiments in that period. Most of the methodology is still valid and should 
be strictly followed, even if this has not always happened in the recent literature. 

The difficulties to correctly interpret the experimental data were partially explained by the observation of Wiegman 
and ter Stege [3] first, and later of Bertotti et. al who noticed as the statistical properties of the noise typically 
vary along the hysteresis loop. Only considering a region where the DW motion is the dominant or better the unique 
magnetization process, as around the coercive field, it is possible to obtain a stationary signal and to correctly perform 
a statistical analysis. This observation t rigg ered a new series of detailed experiments, which were quite well described 
by a new DW model, known as ABBM |l6l IT^. Even if still phenomenological, this model has the merit to describe 
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the statistical properties of the noise on the basis of some simple assumptions verified experimentally, and to be 
analytically tractable. 

This new experimental determination of the noise properties put in evidence a series of interesting characteristics. 
The BK noise is self-similar, and shows scaling invariance and power laws. In other words, it has the typical features 
of a critical phenomenon, whose detailed nature has been the subject of an intense debate. Following this line of 
interpretation, in 1991, Cote and Meisel [islIT^ claimed that the Barkhausen effect is a particular example of a "self- 
organized critical" system '2^, in which there is no need to fine tune a control parameter, as in usual second-order 
transitions. Following this suggestion, new experiments and models were published quite regularly in the literature. 
After nearly 15 years, it is necessary to review this large production of data and ideas and, taking into account the 
important result previously appeared in the literature, to propose a single framework to interpret the phenomenon. 
This is particularly important in order to deal with recent observations of the Barkhausen noise lacking a coherent 
interpretation, such as for measurements in thin films, hard materials, and low dimensional systems. 

II. EXPERIMENTS AND MATERIALS 

In the vast literature about the Barkhausen noise, there is such a large production and analysis of experimental 
data that there is a serious risk of confusion and misunderstanding. Data have been collected on nearly all possible 
soft magnetic materials, and considering many different properties, so that one can, in principle, assume them as 
sufhciently exhaustive and reliable. The absence of a coherent interpretation could then be ascribed to the lack of a 
proper theoretical model. 

Paradoxically, the situation is very different. Despite the large production of data, the number of reliable results 
is quite small, and is restricted, with only a few exceptions, to last decades of studies. There are many reasons for 
this, as we will discuss in detail in the following. One of the main points is that there is no general agreement about 
a 'standard' experimental setup which would allow for a meaningful comparison of the results. In fact, different 
experimental procedures can all produce "nice" Barkhausen data, but their interpretation could widely differ. The 
way the measurement is performed influences significantly the details of the magnetization process. Thus, there has 
been a strong effort in the last years to set experiments in the simplest and most controlled way, avoiding unnecessary 
complications. This opened the way to a convincing theoretical interpretation based on relatively simple models. It 
is only through well defined experiments that parameters can be estimated and used in the theory. 

As discussed in the introduction, a crucial observation is that the noise properties, such as the power spectrum, 
change along the hysteresis loop |l4l |l5j . As a consequence, it is important to measure the noise only around the 
coercive field in the steepest part of the loop, where domain wall motion is the dominant (and often the unique) 
magnetization mechanism, and the Barkhausen signal is stationary. This point questions the reliability of many 
experimental results previously reported in the literature. 

It is convenient to critically discuss the experimental procedures proposed in the recent literature, in order to 
identify a series of practical rules to establish a 'standard' setup. This will enable the direct comparison of different 
experimental data, and the development of a common theoretical framework in which to analyze and understand the 
magnetization process. 

A. Experimental setup 

As shown in the introduction, earlier experiments, starting from the original one by Barkhausen, are all based 
on the induction of a flux change into some concatenated coils in response to a slowly varying external field. The 
regularity of the applied field contrasts with the irregular character of the induced pulses, as a result of the random 
motion of the DW. More recent experiments make use of magneto-optical methods based on the Kerr effect (MOKE) 
and are more suitable to investigate the noise properties of thin films, as we will describe in detail (see Sec. IIII 

1. Inductive measurements 

The essential elements of a inductive Barkhausen noise measurement are simply identified: an external solenoid or 
a Helmholtz coil, able to produce a sufficiently homogeneous field along the sample, and a certain number of pick-up 
(secondary) coils wound around the sample to detect the induced flux. In general, it is a good practise to use a 
solenoid much longer than the sample to ensure a uniform applied fleld. Also in the case of application of an external 
tensile stress, it is preferable to use a non-magnetic ribbon glued to the sample and attached to a weight (as done in 
Ref. pH using amorphous materials), so that the sample remains in the center of the solenoid. To make measurements 
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FIG. 1: Schematic representation of a setup for inductive measurements. In this case, the sample is longer than the solenoid, 
so that the applied field is not uniform. The flux-closing yoke is used to change the apparent permeability [From [l^ . Fig. 1, 
p. 2909] 



at different permeability there are two alternatives with both advantages and disadvantages: as in |22l |. a magnetic 
yoke with a variable gap can be used, but, obviously, with a sample longer than the solenoid (see Fig.^; otherwise, 
the sample can be progressively cut, as in Ref. plj . 

The induced signal <i> both detects the contribution of the applied field H across the pick-up coil (having a cross 
section Apickup) and of the material itself, so that $ = ApickupfJ-oH + AI , where A is the material cross section, and / 
the magnetization. This expression is valid as long as the eddy current shielding is negligible, as happens for the usual 
measurements performed at very low frequency (usually up to no more than a few Hz). In high permeability magnetic 
materials and for A not too small in respect to Apickup, the contribution of the external field can be neglected for 
frequency up to several kHz, thus the induced flux contains all the information about the material properties. On 
the contrary, with low permeability materials, or when A <C Apickup, as for instance, in thin films, it is necessary 
to compensate for the induced flux in air in order to have $ ~ AI. This is usually performed using another pickup 
coil with the same number of turns and the same cross section, but wound in the opposite direction. A more refined 
solution requires a further coil put in air to get a finer compensation 14, 23] . This solution is particularly important 
when it is necessary to measure the hysteresis properties of thin films using inductive methods . 

In the simplest experimental conditions, as for a material having only two opposite domains and thus a single 
DW, the induced fiux is simply proportional to the DW velocity w, as $ = N{2Isdv), where Is is the saturation 
magnetization, d the sample thickness, and N is the number of turns. In the more common case of a multiplicity 
of DWs, one can assume that v measures the velocity of the active DWs, even if this must be considered just as a 
rough approximation. In fact, the induced signal indirectly collects, in a very complicated and unpredictable way, 
the changes coming from all the rest of the sample because of long range spatial interactions, and/or propagation of 
local deformations of DWs. This latter effect has been quantitatively studied, showing that the flux variation decades 
roughly exponentially with a typical scale of the order of a few cm . This suggests to keep the width of the pickup 
coil as small as possible (a few mm, in practice). A measurement obtained with N turns is thus a good approximation 
of a measurement in a single cross-section of the material, with the advantage of the amplification of the signal by the 
factor N. There is also an additional effect suggesting to keep the width of the coils small. The magnetostatic field 
in an open sample produces a counterfield known as demagnetizing field, which essentially depends on the sample 
shape and on the domain structure. This field is spatially constant only in a sample having the form of ellipse and 
at saturation, otherwise it varies, as for instance, when a domain structure is present. The demagnetizing effect has 
a strong effect on the DW dynamics, limiting for instance the maximum size of a jump (see Sec. IIII A^ : it is thus 
convenient to consider a constant demagnetizing field, and this is again obtained limi ting the width of the secondary 
coils. A notable example not following this rule is the setup used by Spasojevic et al. j25| where the pickup coil has a 
length of 5.5 cm, while the sample is only 4 cm long. This seriously questions the interpretation of the data and the 
obtained scaling exponents given the strong spatial variation of the demagnetizing field. 

Even following all these experimental prescriptions, we must clarify that the inductive measurement is always a 
collective detection of many moving DWs, where propagation and long range interactions still play a significant role. 
We emphasize it as the most successful models proposed in the literature (see Sec. IIVII are based on the dynamics of a 
single domain wall. It is thus necessary to clarify what is the reason of their success and, at the same time, understand 
their limitations. In particular, we will describe in detail how the open problems, still under severe theoretical and 
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experimental investigation, could be possibly understood introducing more complicated magnetization models which 
have not been considered so far. 

A strict consequence of this collective detection of DWs motion is that the measured signal cannot distinguish 
between single BK jumps and a superposition in space and/or time of them. This limitation has strong consequences 
on the interpretation of the statistical distributions of BK jumps, and in particular on the effect of the external 
field rate (Sec. IIII Ap . Quite reasonably, it has been noted that at zero driving rate (the so called adiabatic 
limit) the signal should be considered a sequence of Barkhausen avalanches, i.e. of elementary jumps without any 
superposition. On the contrary, at finite driving rates, the signal appears a sequence of pulses containing a certain 
number of avalanches. We adopt this distinction, as it is necessary to understand theoretically the effect of the 
driving rate on the scaling exponents. Clearly, there is a practical difficulty to establish the low frequency regime 
which approximates the adiabatic limit. 

In order to further simplify the measurements, we finally note that the applied field has usually a triangular shape 
instead of the common sinusoidal waveform used in hysteresis loop measurements. As the noise is taken around 
the coercive field, as said, this ensures to have constant magnetization rate, which is the fundamental parameter to 
describe the variation of the BK statistical distributions fSec. and of the power spectra fSec. IIIIB)| . 

2. Magneto-optical Kerr effect measurements 

As described ab ove, earlier measurements on thin films were done using inductive setups, with fine compensation 
of air flux [T^ I23L . These measurements focused both on the noise power spectra , and on the avalanche 
statistical distributions, detected along the hysteresis loop Q^. The authors were able to observe the noise in films 
down to a thickness of 40 nm. This is a very surprising result, as we must considered that the size of avalanches 
roughly diminish with the sample thickness, and thus can be considerably small at such thickness. In ref. [3 the 
largest avalanche size detected is of the same order of magnitude (~ 10~* Wb) found in ribbons (cfr. Fig.^J, while 
the largest duration is two order of magnitude smaller. This means that the 83-17 NiFe film used has DWs with 
an extremely high mobility. In other materials, the situation cannot be so favorable. For instance, we have recently 
used a similar setup to detect the noise in Co-base thin films. For thickness below 1 /im, the BK signal becomes 
undistinguishable from the instrumental background noise, so that the measurement is practically impossible. Thus 
the only valid alternative is to use optical methods, based on the magneto-optical Kerr effect (MOKE), which is the 
change of polarization of the light reflected by a magnetic material. This technique is extensively used to investigate 
the domain structure of most magnetic materials |2q . There are fundamentally two variations of the basic MOKE 
setup applied to the investigation of the Barkhausen noise and in particular to the statistical distributions. The first 
was introduced by Puppin et al. who added an optical stage to tune the laser spot size onto the film, with 

values ranging from 20 fim up to 700 ^m. The reduction of the background noise is obtained using a photoelastic 
birefringence modulator working at 50 kHz, while the signal is detected with a lock-in amplifier at 100 kHz to enhance 
the signal-to-noise ratio. This system allowed to estimate the distribution of avalanche sizes in a Fe thin films IsollsH 

(see sec inrnii . 

The second solution makes use of an advanced video processing technique to directly observe the domains and 
their motion Is^ . Is^ . The main instruments used are an objective of numerical aperture of 0.9 having a xlOOO 
magnification. The domain images are obtained at the rate of 30 frames/s in real time and a spatial resolution of 400 
nm |33j| . In this particular setup, the external field is applied with an electromagnet up to nearly 99% of the coercive 
field starting from the saturation. Barkhausen jumps are thus triggered only by thermal fluctuations, so that any 
sweep rate effect is completely eliminated. The experiment is performed repeatedly imaging the same film area, using 
up to 1000 repetitions. The BK jump size is defined as the area having two end points outside the image; this choice 
raises some problems in the interpretation of the experimental data, as we will discuss later fSec. IIIl'CI) . 

B. Magnetic materials 

Nearly all soft magnetic materials have been studied using the Barkhausen noise technique. But despite the large 
amount of published papers, a systematic identification of material properties appears intrinsically hard to perform, 
because of the large differences between the same type of material, the possible annealing procedures, etc. Up to 
now, the only identification of different classes of behavior (the "universality classes" deeply discussed in Sec. IIV D 2|) 
has been reported for SiFe polycrystalline ribbons with high Si content and amorphous FeCoBSi alloys j2ll |. which 
show two distinctive scaling exponents ruling the Barkhausen distributions. This difference is understood in terms of 
a different DW dynamics, as it will be discussed extensively. 



L : stabilized He-Ne laser 
P : polarizers 
B : beam expander 
F : best form focussing lens 
C : coils for magnetic field generation 
vS : sample 
FL : focussing lens 
M ; pholoclastic modulator 
A : second polarizer 
D : silicon photodiode 

FIG. 2: The focused magneto-optical Kerr magnetometer with variable laser spot sizes used to investigate the avalanche size 
distributions in thin films [From [S^, Fig. 1, p. 1753] 

Earlier fundamental studies have been performed in frame- type sing le crystals to eliminate the effect of the de- 
magnetizing field, where a single DW is present (see Ref. 0, HE l^raTIs^ for 3% wt. SiFe, ref. 4 for 3.8% SiFe, 
Ref. i ^ for 4% SiFe, and Ref. 10 for 4.5% and 6.5 % SiFe). The most studied polycrystalline materials are the 
traditional non-oriented 3% wt. SiFe strips [H [Hill EH , or the 1.8% ones as well as the GO SiFe alloys 
[Til HE 13, liE EE E3, El • other interesting SiFe alloys are the ribbons (thickness about 50-60 /im ) with high Si 
content (% 6.5-7) produced by planar flow casting, which reveal interesting noise characteristics .21. 49, 50.1. A large 
attention has received also the NiFe alloys, such as the Permalloy (Ni8oFe2o) and similar 0, Il4l 1411151115^ or the 
Perminvar (30 % Fe, 45 % Ni, 25 % Co) |53l l54l|. particularly when a strong disorder is induced by precipitation of 
some crystalline phases. 

In more recent years, many studies have been focused on amorphous materials, especially the iron-based family: 
some commercial ribbons, like the Metglass 2605 El El EE Ei, the Metglass 2826 or the Vitrovac 6025 Igi, and 
the vast production of various laboratories with different composition and additiv es Iplllisl 11^ Isollssl 15^16011611 . On 

the contrary, Co-based amorphous alloys have received a very limited attention j4ll |. due to their high permeability 
values connected to the usual presence of a single large BK jump. 

Despite the large variety of compositions used, we should stress here that not all the soft magnetic materials are 
appropriate for BK measurements. Especially, materials having very nice 'soft' properties (high permeability, small 
losses, etc.) show reduced or no BK noise, like the Co-base amorphous alloys or a well prepared Permalloy, with 
limited structural defects. In fact, the essential ingredients for a 'nice' noise are the presence of a consistent structural 
disorder, hindering the motion of DW, and a non-negligible demagnetizing effect, limiting the size of avalanches. Also 
the thickness appears to be an important factor to take into account, as it sets the limits of eddy current shielding 
effect. For instance, SiFe ribbons having thickness of the order of hundreds of /im, have to be magnetized at very low 
frequencies to detect well defined BK jumps [l^lioj . 

III. EXPERIMENTAL RESULTS 
A. Barkhausen signal and avalanches distributions 
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In Fig. 121 we show some typical experimental BK data obtained in a polycrystalline Fe-Si 7.8 wt.% ribbon and in an 
amorphous Fe64Co2iBi5 ribbon, the latter measured under moderate tensile stress. The different sequences, acquired 
during the same time interval, correspond to increasing applied field rates. At low rates, the BK signal clearly comes 
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FIG. 3: Time sequences of Barkhausen noise measured in a polycrystalline FeSi 7.8 % wt. ribbon (left) and an amorphous 
Fe2iCo64Bi5 under moderate tensile stress (right). The labels indicate the frequency of the applied field. 

out as a sequence of distinguishable and well separated avalanches: at larger rates, this separation is progressively 
lost and the noise resembles a continuous sequence of peaks. This feature is common to both materials, even if this 
similarity is only seeming, as their statistical properties are very different and appear to be representative of two types 
of behavior, the 'universality classes', with different DW dynamics, as we will describe in detail. 

To investigate the statistical distribution of BK jumps is necessary to univocally define the size and the duration 
of an individual jump. This poses a delicate problem, given the self-similar fractal nature of the BK signal, and the 
fact that, in addition, the background noise limits the detection of smaller avalanches. The common way to solve 
this problem is to define a threshold, a 'resolution coefficient' 0, H^, |^ 113, E3 1 in order to set the temporal hmits 
of a single avalanche. This procedure has been justified in the context of the fractal properties of the signal, as it is 
possible to associate to it a random Cantor dust, with a fractal dimension < 1, which can be used to calculate the 
critical exponents . The determination of the threshold can vary, but in general we found the results are nearly 
insensitive to changes in its value |6^, as long as it is not too small or too large. Because the BK signal amplitude v 
follows a power law of the type 



Piv) 



-(l-c) 



fiv/vo) 



(1) 



where c is proportional to the driving rate, we usually assume the threshold between 5 and 15 % of wo. 

With this definition, the duration T and the size 5" (the integral within the time limits) are unambiguously defined. 
Like the signal amplitude, also the distributions of T and S follow power laws 



p(T) = r-"g(r/To), 



(2) 



and 



P{S) = S-^fiS/So), 



(3) 
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Ref. 


Material 


Type 


Dimensions 


r 


a 


1/avz 


Rate dep. 


m 


81%NiFe 


Wire 


50cm X <1> 1 mm 


1.73 


2.28 


1.63 


Yes 




Vitrovax 6025X 


Ribbon 


4 cm X 1 cm x 30 /im 


1.77 


2.22 


1.51 


? 


m 


Perminvar 


Strip 


5 cm X 2 cm x 100 nm 


1.33 






No 


Talk 


Annealed steel 


Strip 


2.5 cm X Ij cm x 840 


1.24;1.27 






iNo 


@ 


Fe64Co2iBi5, 
Fe2iCo64Bi5 
amorph. (tensile 
stress) 


Ribbon 


28 cm X 1 cm x 20 /im 


1.3 


1.5 


~1.77 


No 




Fe2iCo64Bi5 
amorph. (no stress) 


Ribbon 


21 cm X 1 cm x 20 ^m 


1.46 


1.74 


1.70 


No 


m 


SiFe 1.8% 


Strip 


20 cm X 1 cm x 180 /im 


1.5 


2 




Yes 


m 


SiFe 6.5%, SiFe 
7.8%, Fe64Co2iBi5 
partially cryst. 


Ribbon 


28 cm X 0.5 cm x 48 fim 30 cm 
X 0.5 cm X 60 fim 30 cm x 1 
cm X 20 j-im 


1.5 


2 


~2 


Yes 



TABLE I: Experimental critical exponents reported in the literature which can be considered sufficiently accurate and reliable. 
Exponents r, and a are the critical exponent of the size and duration distributions, respectively. The exponent 1/auz relates 
the average size of an avalanche to its duration (see text). "Rate dep." indicates if the exponents r and a depend on the applied 
field rate. 

where Tq and 5*0 are two cutoffs, and a and r the critical exponents. 

It is interesting to note that, despite the large number of papers regarding the BK noise, the reliable estimates 
of these critical exponents are very limited. We report in Tab. HI the results presented in the literature we believe 
are sufficiently accurate to be considered and analyzed. We leave out, for instance, the results of the quite famous 
papers of Cote and Meisel Eliaj, 

who first suggested the BK exhibits self-organized criticality 20], because they 
are based on a very poor statistics. At first sight, the measured exponents span a quite large range of values, and no 
general behavior can be easily guessed. But it is worth to consider in detail each experiment analyzing the particular 
conditions of the material and of the measurement to find a possible explanation for this large variability. In particular, 
two key factors appear to affect the results, and should thus be carefully checked: the applied field rate dependence 
of the critical exponents, and if the data are taken in a small bin around the coercive field or along all the hysteresis 
loop. After this warning, we proceed to analyze the literature results in more detail: 

1. The 81%NiFe used in Ref. :51 is the only wire which has been measured accurately. The authors use 3 different 
samples with different annealing conditions, but do not report any indication about the structural characteristic 
of the materials. They also did notice a dependence of the statistics on the applied field rate, but to simplify 
matter, they made the measurements at the same magnetization rate, which later has been established to be the 
correct parameter to take into account |l6l Il7l | . As it is not explicitly mentioned, the noise has been probably 
recorded along all the hysteresis loop and it is thus likely that it is not statistically stationary. 

2. Similar critical exponents have been reported by Spasojevic et al. p5l |. with a completely different material. As 
mentioned in Sec. Ill Al this experiment has some serious drawbacks due to the large extension of the pickup coil 
in respect the sample length, which have been shown to affect the size distributions @ . The data are taken in a 
small field interval, presumably around the coercive field, but it is not clear from the paper if the applied field 
is large enough to saturate the material. The rate dependence of the critical exponents has not been discussed. 

3. The Perminvar used in (H^ has been subjected to an annealing at 450 °C for 24 h. This prolonged treatment has 
the effect to induce the precipitation of micro-crystalline phases, which create a large number of strong pinning 
centers for the DW, resulting in a low permeability of about 250. The estimation of the critical exponent is 
limited to the size exponent t, but it is made with high accuracy. 

4. A comparable size critical exponent has been measured for a sheet of interstitial-free, low-carbon, annealed steel 
j64j. The external field is applied both parallel and perpendicular to the rolling direction. The corresponding 
hysteresis loops show a different permeability, which affects the value of the cutoff 5*0, and slightly the critical 
exponent r. As matter of fact, the estimation at the lowest permeability (smaller Sq) appear much less accurate, 
due to the limited region of the power law dependence. 

5. Ref. 43] reports the first estimation of the critical exponents in a polycrystalline SiFe strips (1.8 % wt.). Within 
the experimental errors, the critical exponents show a linear dependence on the applied field rate, expressed by 
a = 2 — c, and r = 1.5 — c/2, where c is the parameter comparing in Eq. ^ 
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FIG. 4: (Left) Distributions of Barkhausen jump sizes measured in different materials. The solid line has a slope r ^ 1.5 while 
for the dashed one r ~ 1.27. (Right) Similar plot for duration distributions. The solid line has a slope a ~ 2, while for the 
dashed one a ~ 1.5. [From Fig. 1, p. 4706] 

6. In Ref. l2lL we have reported measurements of the critical exponents in polycrystaUine SiFe ribbons (with high 
Si contents) and amorphous alloys with composition FcxCoss-xBis. We also measured a partially crystallized 
Fe64Co2iBi5, obtained after annealing for 30 min at 350 °C and then for 4h at 300 °C under an applied tensile 
stress of 500 MPa. This induces the formation of a-Fe crystals of about 50 nm, with a crystal fraction of ^ 5% 
[67l |. The SiFe alloys and the partially crystallized alloy show r ~ 1.5 and a ^ 2, measured at the lowest 
possible frequency / — 3-5 mHz (see Fig. 01- left). They also show a linear dependence on the applied field rate, 
as shown in Fig. [SJleft ^3)E3- The highly magnetostrictive Fe^^Coss-^Bis alloys (A^ ~ 30 — 50 x 10^^) have 
been measured under a tensile stress of cr '--^ 100 MPa. The applied stress is found to enhance the signal-noise 
ratio, reducing biases in the distributions, but does not change the exponents j68l |. as long a small stress is 
applied (see below). These alloys yield r ^ 1.27 and a ^ 1.5 (Fig. ^ right), independent of the apphed field 
rate, as shown in Fig. [SJ-right. In contrast, the shape of the cutoff changes drastically as the rate increases. 

7. The authors of Ref. '6^ used the same type of amorphous alloy of Ref. l2ll but without the application of a 
tensile stress. The resulting critical exponents are larger, and the authors have no a simple explanation of 
the difference. As they claimed, the residual stress may be the origin of the anomaly. We can confirm that 
unstressed amorphous materials have in general hysteresis loops with restricted regions of linear permeability, 
thus the determination of the exponents can be difficult and less accurate. 

In Tab. ^ we also report the values of another critical exponent \/(jvz relating the average size {S) to the duration 
T as (SV- yi/^i^. This exponent has been shown to be in strict relation with the power spectrum exponent, (see 
Sees. HiTbI and rvUl. 

To understand the origin of the variability of the critical exponents is important to first focus on the values which 
has been found to be close in different samples and experiments. In addition, it is better to consider the data recorded 
in a small bin around the coercive field, that can thus be considered as stationary. Two sets of exponents emerges: 
one is measured in polycrystaUine alloys or with embedded grains, characterized by t « 1.5, and a w 2, with a 
linear dependence on the applied field rate. The other is measured in materials where the crystalline matrix is absent 
(amorphous samples under tensile stress) or progressively lost by induction of a strong disorder (like in the Perminvar 
of Ref. 113). They show r ~ 1.3, and a w 1.5, and no rate dependence is found. These results support the idea 
that there are (at least) two kinds of behavior, the universality classes, to which correspond different sets of critical 
exponents. We will came back to this matter in more detail in Sec. IV Al where we will also try to explain the 
occurrence of the other experimental exponents. 

Given the variability of the number shown above, it is quite clear that the sole comparison of the values of the 
exponents r and a can be misleading. In other words, it necessary to determine other experimental critical exponents 
which can confirm the existence of the two universality class. In this respect, an important suggestion has been 
given by Perkovic et al., who pointed out that "an explanation for the experiment must involve collective motion of 
many domains; it must provide an explanation for the power-law scaling regions, an it must provide an explanation 
for the cutoff [From 69], pg. 4528]. We have already stressed the importance of the demagnetizing effect in the 
DW dynamics, so that we investigated in detail the effect of the change of the demagnetizing factor k on the power- 
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FIG. 5: Distributions of the avalanche sizes S for the polycrystaUine SiFe (left) and the amorphous (right) material at different 
applied field frequencies. In the polycrystaUine sample, the critical exponent r varies linearly with the frequency [43 . |6^ 
between 1.5 and 1.15, while it is independent in the amorphous one r ~ 1.27. 




FIG. 6: The cutoff of the Barkhausen size (a) and duration (b) distributions as a function of the demagnetizing factor k in 
Fe-Si 6.5 wt% alloy (empty symbols) and an amorphous Fe2iCo64Bi5 sample under constant tensile stress (filled symbols). 



law distributions. We used two samples belonging to the different universality classes, and measured the avalanche 
distributions, progressively cut the samples to change the apparent permeability ^app, and thus calculate the factor 
k using the well known expression fj^app = fJ-fi^ + k^), with /x the intrinsic permeability. We verified that the cutoffs 
So and Tq strongly depend on k, as in Fig. (HJ and follow the scaling relations Sq k~^^'^'', and Tq ^ k~^*'. Best fit 
estimation of the exponents gives 1/ak ~ 0.57, Afc ~ 0.30 for the polycrystaUine class, and l/ak ^ 0.79, Afc ^ 0.46 
for the amorphous one. 



B. Power spectra and avalanche shapes 



If we merely count the huge number of papers in the literature concerning the spectral properties of the BK noise, 
the contrast with those considering scaling exponents results clearly. Since the earlier papers, but especially during the 
70s-80s, almost all the authors publishing about the BK effect considered the problem of the detection and explanation 
of the power spectrum. The main reason for this abundant production is related to the practical need to reduce the 
spectral noise in applications where soft materials are used. In addition, it seemed easy to explain all the details 
the power spectral shape in terms of a superposition of 'elementary' Barkhausen jumps (see for instance Q). These 
elementary units were generally defined a priori, for instance, as exponential or square events: a proper superposition of 
these independent jumps can give practically any measured power spectrum shape. Even if successful to give a general 
description of the power spectrum shape, this approach was not able to account for the microscopic magnetization 
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events and the details of the DW dynamics, so it really did not represent a significant step towards the fundamental 
understanding of the BK noise. Remarkably, after the publication of the two seminal papers by Alessandro et al. 
0,13 1 introducing the ABBM model fSec. lIVB)) . the investigation of the power spectrum apparently lost interest, 
even though the problem remained basically a puzzle. Only recently the subject, seen under a completely new point 
of view, received a new impulse. 

As it should be clear from the discussion in the preceding sections, the detection of stationary or non-stationary 
Barkhausen signals has a strong effect also on the estimation of spectral properties. More precisely, traditional power 
spectral calculations using FFT must be applied only to stationary data, otherwise they can give unpredictable results. 
A non-constant value of the permeability, for instance, introduces low frequency components unrelated to the BK 
noise. Thus we must consider only two types of spectral estimation: i) those taken with a single DW at constant 
velocity, as in [Tolls^ IstL l70l | . where a feedback setup keeps the flux rate constant, and ii) those taken in a small bin 
around the coercive field, as first performed in 

Measurements of power spectrum of type (i) are particularly important. In fact, not only they are a measure of 
the time correlation of the BK signal (in the frequency domain), but more they are a measurement of the spatial 
fluctuation of local fields which are needed to unpin the DW. In other words, these measurements describe the spatial 
spectrum of the random pinning field (see Sec. IIVB|I . Measurements on single crystal of SiFe give a power 
spectrum, and thus a k^^ spectrum of the pinning field [lol l36j| . because of the linear relation between / and k at 
constant v. It is curious to note that measurements on a toroidal Metglass sample gave a k~^-'^ spectrum Even 
if a single measurement cannot be taken as conclusive, this result can be regarded as a manifestation of a different 
universality class, similar to what reported for the avalanche critical exponents. 

As anticipated, the most accurate report of power spectra in a polycrystalline material has been presented in 
Ref. A strip of non-oriented 3% SiFe, having a grain size of about 100 fiui is measured under a small tensile 

stress {'^ 5 MPa), and at variable apparent permeability obtained varying the air gap with a mumetal yoke (see 
Fig.0. The effect of the grain size on similar materials is reported in Ref. [jy. Successively, a few other papers have 
reported the power spectra of single crystal and polycrystalline SiFe materials [l^l , as well as of amorphous materials 

There some general common characteristics that can be summarized as follows: 

1. The spectrum F{uj) has a typical shape at high frequency, with — 1.7 ^ 2, and scales linearly with the 
average magnetization rate I, so that the spectra normalized with the average flux F{u))/SI coincide at high 
frequency (Fig. [TJ; 

2. At lower frequency, the spectrum shows a marked peek at a frequency roughly proportional to /^/^, with an 
amplitude which takes a constant value at low rates, as shown in Fig.[7| This amplitude roughly scales with the 
permeability, as /i^'^ W^ - 

3. At frequencies lower than the peek, the spectrum scales as , with -0 ^ 0.6, or ~ 1. 

These general features describe only qualitatively the observed results. More precisely, there are at least four 
different typical shapes, as shown in Figs. 18191 

Single crystal and polycrystalline SiFe strips show the simplest power spectrum, with 'd 2, and ■0 0.6. The 
polycrystalline SiFe ribbons at high Si content keep the same exponents, but present a large region at intermediate 
frequencies where the spectrum strongly deviates from the simple 1//^ law. Amorphous materials under tensile 
stress show instead d — 1.7 1.8, and ■0 ~ !> as shown in Fig. O Similar data are found in the unstressed case, 
■d — 1.73 ± 0.08 and i/' ~ 1, for the same amorphous alloy 66] (private communication). Amorphous materials, where 
a proper annealing induces a partial crystallization, show a sort of intermediate case: the high frequency part is more 
similar to the polycrystalline SiFe samples with § 2, while ■0 ~ 1 keeps the values of a pure amorphous material. 
The former property is not a merely coincidence, as we observed a corresponding change in the critical exponents of 
size and duration distributions towards the values measured in SiFe, as shown in Tab. ^ Curiously, all the spectra, 
despite these differences, show the same dependence of the peek frequency on /^/^. 

The interpretation of all these complicated details still represents a challenging puzzle, although some progress 
has been recently made |6^[3- Earlier estimates of the exponent 'd (HJ, and similarly in jl^E^jIzll) led to the 
conclusion that "the asymptotic decay of the spectrum is not given by the shape of the pulses, but is exclusively 
determined by the distribution of the pulse duration and the correlation between the size and the duration" [Adapted 
from 51], pg. 15]. In our notation, '0 = 1 + 2avz — a, or analogously (3 — T)/avz (see also Sec. lVB|l . Using the data 
of Tab. ^ this implies i9 ~ 3, much larger than the experiments. This result was successively found to be valid only 
for r > 2, while for r < 2, Kuntz and Sethna [t^] calculated i? = 1/avz, which is the same critical exponent relating 
the average avalanche size to the duration (Tab. P). This prediction agrees pretty well with the experiments, as we 
can see in Fig. 1101 In particular, it works for the amorphous material in an extended frequency range; the same is 
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FIG. 7: (Left) Normalized power spectra F{uj)/SI at fixed permeability of a poly crystalline 3% SiFe strip. The solid lines are 
the experimental results while the dotted lines are the prediction of the ABBM model fSec. lf V HI : (right) Normalized maximum 
of the power spectra Fm/SI for permeability ranging from 2500 (o) up to 60,000 (t>). Continuous lines are guide to eyes only. 
The bold line is the boundary between the intermittent (c < 1) and the continuous (c > 1) DW motion. Here c is equivalent 
of the exponent c of eq.[T] [From Ref. [l^], figs. 7 and 10, pg. 2912-2913] 



not valid for the polycrystalline SiFe ribbon, and the reason will be discussed after introducing the hypothesis at the 
base of the calculations by Kuntz and Sethna. 



1. Avalanche shapes 

The simple result — \juvz obtained by Kuntz and Sethna [T^I is based on the existence of a certain number 
of scaling relations and on a few hypothesis, like the complete separation of the avalanches in time. We leave the 
discussion of these important aspects to a further section (see Sec. IVB|I . focusing here on the new proposed scaling 
relations which can be tested experimentally. These relations involve universal scaling functions whose shape could 
provide a critical test for model, as "a sharper tool for discriminating between different universality classes than 
critical exponents" . 

The first relation states that the average avalanche shape scales in a universal way so that 

vit,T)=T'/''-^-'Uapeit/T), (4) 

where v is the BK signal, t is the time and fshape{t/T) is the universal scaling function. Similar to this relation, we 
can consider the voltage as a function of magnetization s — vdt to get 

v{s,S) = S^-'''''gshape{s/S) (5) 

A third relation deals with the fluctuations of avalanche sizes, and considers the probability P{v\S) that a voltage 
V occurs in an avalanche of size 5*. This probability scales as: 

P{v\S) = V'^f^oltageivS^-''-^) (6) 

and, again, f voltage is another universal scaling function. 

A partial experimental verification of these relations has been reported in a couple of papers [65l l66j . In Fig. 1111 we 
show the complete universal scaling laws of Eqs.2|and|Slfor the same 7.8 % SiFe ribbon and the amorphous sample 
of Fig. 13 scaled using the theoretical estimation of the critical exponent avz. The experimental data are compared 
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FIG. 8: Normalized power spectra of a single crystal 3% SiFe strip (left, reproduced from [t^ . Fig. 3, pg. 580], and a 
polycrystalline 7.8 % SiFe ribbon (right) as a function of the magnetization rate I. 




Frequenty (Hzt Frequency (Hz) 

FIG. 9: Normalized power spectra of the same amorphous material of Fig.|3(left), and of a partially crystallized (16%) sample 
induced by annealing (right, reproduced from fsoll . Fig. 2, pg. 300) as a function of the magnetization rate /. In the left Fig., 
/* is roughly the frequency where the shape changes from to 



with the theoretical average shape predicted using the ABBM model (see Eas. llUlllU2ll . In the case of the amorphous 
material, the scaling is pretty good, except for the magnetization signal at small sizes. In particular, both time signals 
show a marked time asymmetry, similar to what found in 25] where the average of avalanches is made without any 
scaling. The asymmetry is also visible in the magnetization signal, even if less marked. Here, the shape calculated 
using the ABBM (see Eg. 11021 and Sec. IV C|) reproduces pretty well the experiments, except for the small sizes in the 
amorphous material. It is interesting to note that all the models presented in the literature predict symmetric time 
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FIG. 10: Comparison of the power spectrum F{uj) with the average size (S) as a function of the inverse of duration T, for the 
amorphous material and the polycrystalline SiFe ribbon (7.8 %, inset) of Fig. 0] The theoretical prediction of Ref. [74| is also 
shown, with 1/aiyz equal to 1.77 and 2 [From |65|. fig. 1, pg. 1086] 



and magnetization signals. In Sec. IV CI we will discuss the relevance of the shape of the scaling function as indicator 
of universality class. 

The scaling of time signal data of 7.8 % SiFe ribbon using the theoretical exponent = 2 is clearly incorrect. 

The reason is illustrated in Fig. where we plot the average avalanche size (S) as a function of duration T. The 
exponent l/avz (the slope of the curve in the log-log plot) slightly varies from ~ 2 down to about 1.75. If we use this 
'variable' exponent l/avz^, the rescaling of the time signal is highly satisfying. 

The third universal scaling function (Eq. Ej) has been reported in j65l | for an amorphous material, as shown in 
Fig. El This confirms the validity of scaling relations proposed in [t^I, at least for the amorphous material. In fact, 
the relation in Eq. Eldos not hold so well for the 7.8 % SiFe ribbon. 



2. High order power spectra 

In Ref. t55i|, O'Brien and Weissman introduced a further statistical characterization of BK signal calculating power 
spectra of order higher than those shown in the previous section. This analysis is motivated by the possibility to prove 
the BK as a real example of SOC '20|, lacking of any inertial effect unlike, for instance, grain falling in sandpiles. 
In fact, the mere "existence of pulses with a range of sizes and durations should not be taken as conclusive evidence 
for SOC" [From |^, pg. 3446]. Thus it is important to understand if the observed scaling properties reflect a real 
self-organization involving "cooperative avalanche processes invoked by SOC", or simply reflect the scaling properties 
of the pinning field, or, alternatively, the vague proximity of a disorder-driven critical transition |69l ]. 

The high order spectra involve correlations between different temporal scales, so for instance, could clarify if events 
on a certain scale systematically precede event on another scale, as in the case of systems showing precursors such 
as earthquakes. Another important aspect, more specific of BK noise, is the possibility to distinguish whether the 
high frequency spectrum primarily comes from short independent pulses, as often assumed, or conversely from 
the fine structure of longer pulses 73]. This is addressed considering the 1.5 power spectrum S'i.5(/2, /i), involving a 
third moment of time signal v{t), defined as: 

^ F{vit)}F{Hit,h)r 



where H(t,fi) is the time-dependent Haar power at the frequency /i, F{-} is the usual Fourier transform those 
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FIG. 11: Average avalanche shapes for the 7.8 % SiFe ribbon (left) and the amorphous sample (right) of Fig. |21 Time signal 
v{t, T) and magnetization signal v{a, S) are scaled according to Eas. l4l5l Bold lines are the theoretical predictions of Egs. llUll 
11021 a semicircle and an arc of sinusoid, respectively (see Sec. IV (Jll . Numbers in the graphs denote avalanche size (in nWb, 
upper figures) and duration (in milliseconds, lower figures). 




duration r(ms) t/7 



FIG. 12: Average avalanche size as a function of the duration for the 7.8 % SiFe ribbon (left). Time signals are rescaled as in 
Eq.|l| using a variable exponent 1/auz* , deduced by the plot on the left. 



components are calculated at frequency /2 |73|. Similarly, the second spectrum S'2(/2,/i) is defined as: 

,,_ F{H{t,h)}F{H{t,h)Y 
S2U2J1)- <H{tJ,)>^ 

A complete presentation of these spectra have been reported in Ref. (jsl for a Fe2iCo64Bi5 amorphous ribbon (in 
the unstressed state). In Fig. ^^we show the 1.5 power spectrum. The weak independence of Re{Si.^{f2, /i)} on /i 
indicates that the high frequency spectrum 1//'' contains a large contribution from pulses with duration much longer 
than I//1. This definitely demonstrates that, in this case, the power spectrum does not come from independent 
short pulses. The imaginary part, instead, is sensitive to the asymmetry of the time signal. The positive sign of 
/m{S'i 5(/2, /i)}/i?e{5i.5(/2, /i)} suggests that, on average, high frequency components precede the low frequency 
ones. This is in accord with the time asymmetry of the average avalanche shape of Fig. 1111 
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FIG. 13: Distribution of time signal at fixed size as predicted by Eq.^lfor the amorphous material using the theoretical scaling 
exponent 1/avz = 1.77. [From Fig. 3, pg. 1087] 
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FIG. 14: The 1.5 power spectrum defined in Eq. Qfor a Fe2iCo64Bi5 amorphous ribbon in the unstressed state, at 5 different 
frequencies /i. [From fT^ . figs. 6 and 9, pg. 6366-6367] 



Interestingly, high order power spectra of 1.8 % SiFe single crystals have a completely different behavior [73 j as 
shown in Fig. Here S'2(/2, /i) is nearly flat as a function of /2 and falls off as fi^'^, which indicates that most of 
the power at /i comes from separated pulses of duration l//i. The 1.5 power spectrum confirms this hypothesis, as it 
strongly depends on /i, in contrast with the weak dependence in the amorphous ribbon. In this case, it is appropriate 
to assume the power spectrum as a superposition of individual independent pulses. 
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FIG. 15: Second spectra S2 and 1.5 power spectra S1.5 for a 1.8 % SiFe sing le crystal [From "t?], Fig. 3, pg. 1172]. 

C. Thin films 

All the results presented so far are related to materials which are essentially three-dimensional, with the sole 
exception of the wire reported in Ref. |^ (see Tab. Q}. The behavior of the BK noise in 2D thin films is much less 
known, and the number of papers reporting a complete set of experimental results is very limited 0, |23, 113, 113 ■ 
Many questions remain to be solved, both on the experimental and the theoretical side, so that the topic represents 
an interesting and promising field of study for the near future. 

Contrary to bulk three dimensional systems with a relatively simple magnetic structure of nearly parallel DWs, 
thin films show richer and often more complicated DW patterns. The film thickness plays a fundamental role because 
of the increasing importance of stray fields in the direction perpendicular to the sample, as first pointed out by Neel 
[78.] . In addition to usual Bloch DWs, other different types exist, such as symmetric and non-symmetric Neel walls, 
and the cross-tie walls, a complicated pattern of 90° Neel walls. A complete description of the existing configurations 
together with many experimental images can be found in the excellent book of Hubert and Schafer '2811 . To have an 
idea of the complexity, it is worth reporting their very words about the role of film thickness: "The typical DW in very 
thin films, i.e. the symmetric Neel wall, needs basically a one-dimensional description of its magnetization structure. 
The stray field of this wall reaches out into the space above and below the film, however, and cannot be treated 
one-dimensionally. Cross-tie walls in thin films need at least a two-dimensional description for the magnetization and 
a three-dimensional treatment of the stray field. Bloch walls in thick films and their relatives are two-dimensional 
both in the magnetization and in the stray field. Most complex are cross-tie walls in thicker films ( 100 nm for 
Permalloy) , which need a three-dimensional description of both magnetization and field and which have never been 
analysed theoretically up to now" [From Ref. ^2§\, pg. 241]. 

It is clear that the experimental BK data can be understood only through a precise knowledge of the active DWs, 
and their role in the magnetization process. Considering that the present research is still at a preliminary stage, 
the task will surely require some time to be fulfilled. Nevertheless, it is worth to discuss in detail the main results 
presented in the literature. A quite impressing work has been done in the Ph.D thesis of N.J. Wiegman in 1979 2^, 
who explored the properties of Permalloy thin films using a inductive setup. A selected number of results had been 
previously reported in a couple of papers 0, . 

The simple visual comparison of a typical noise data of Permalloy films with those of SiFe or amorphous samples 
(see Fig.O reveals striking differences. As shown in Fig.^| the BK avalanches appear as sharp peaks, well separated 
in time, with nearly no inner structure. This behavior is confirmed by the data of the size and duration distributions, 
where the maximum duration is nearly two orders of magnitude less than the corresponding value in ribbons, as noted 
before, while the same does not occur for the maximum size. In other words, the avalanches in these films are rather 
sharp peaks of short duration, reversing a magnetization volume comparable to the value observed in ribbons. The 
behavior of the critical exponents r (ai in Wiegman's notation) and a {61), shown in Fig. 1171 as a function of the film 
thickness is quite intriguing. The exponents do not show any strong thickness dependence, except for the increasing 
fluctuations below 100 nm, which, by the way, is the region where complex cross-tie walls are observed. In the same 
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FIG. 16: A typical time signal for a permalloy thin film [From '2^, Fig. 3.2, pg. 44]. 
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FIG. 17: Values of the critical exponents r (ai in Wiegman's notation, left) and a {Si, right) for permalloy films as a function 
of the thickness. [From [H, Figs. 3.10 and 3.14, pg. 57;62] 



region, the critical exponent l/aiyz (x) displays a sharp transition of from 2 to 1.5, as clearly shown in Fig. 1181 

These interesting results should be considered as the basis for further investigations. In fact, the critical exponents 
have the same values observed in polycrystalline ribbons (t ~ 1.5, a ^ 2, l/ai^z ~ 2) for thicknesses down to 100 
nm, suggesting a typical three-dimensional behavior with predominant long-range interactions. The situation at 
smaller thicknesses is more complicated and the role of DW pattern, such as cross-tie walls, is not clear. While an 
accurate estimate of the critical exponents is desirable, it is also necessary to perform an extended scaling analysis. 
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FIG. 18: The critical exponent l/ai/z (x in Wiegman's notation) for permalloy films as a function of the thickness. [From |2^ . 
Fig. 3.7, pg. 53]. 
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FIG. 19: Distribution of magnetization steps AAf in a 90 nm Fe film. Data taken at different laser spots are properly rescaled 
to show a single distribution [From ,30^1 . Fig. 5, pg. 5418]. 
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FIG. 20: (a) Domain images of same 400 x 320 fj,m^ area for a 25-nm Co film between and 4 sec. The sample was 
saturated downward first, and then a constant field was applied upward, denoted by the solid arrow, during observation, (b) 
Magnetization reversal curves obtained from the corresponding domain patterns of (a) [From 32], Fig. 2, pg. 087203-2]. 



not restricted to a single exponent (i.e. t as in most recent studies). 

In Ref. [s^l, the critical exponent t is estimated in a Fe film of 90 nm using the MOKE setup (Sec. Ill A|l . and 
considering different laser spot sizes from 20 to 700 /zm. A single distribution P{AM) is obtained properly rescaling 
the size of the avalanche AM and the distribution P, as shown in Fig.^| Unfortunately, data at small spot sizes are 
not averaged over different positions, so that the overall distribution appears quite unreliable under 100 /im. Given 
that, the estimated exponent is r = 1.1 ± 0.05. A successive estimation using a spot of 20 /im averaged over 10 
different locations yields r = 1.14 [sj. It worth noting that these values are much smaller than the ones observed in 
permalloy, a fact we will comment very soon. The other critical exponents are not considered, even if their estimation 
appears possible. 

The paper of Kim et al. reports an estimate of the critical exponent r in Co films, using a different MOKE 
setup able to visualize directly the DWs. As shown in Fig. 1201 both ends of the DWs are usually outside the frame, so 
that the authors deliberately neglect those jumps having one of the ends inside the frame. The resulting power-law 
distribution has an exponent close to 1.33, averaged over 4 different thicknesses, from 5 to 50 /zm. As we explain 
more extensively in Sec. IV Dl this exponent is different from r, which should be evaluated, in principle, considering 
only the jumps fully inside the frame. As a matter of fact, this appears quite difficult to perform, given the limited 
number of such jumps. 

It worth noting that a similar comment could be appropriate also for the previous experiment of Puppin. In that 
case, it is not possible to distinguish between jumps fully inside the laser spot and the others. Moreover, the relative 
probability of the two types of events is not known, so that in principle we cannot be sure if the exponent t is correctly 
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estimated. With many others, this could be one of the reasons to get the lower value in respect to the permalloy. 
But, as discussed above, all these matters have be confirmed by further experiments. 

IV. MODELS AND THEORIES OF MAGNETIZATION DYNAMICS 

A theory of the Barkhausen effect should explain the statistical features of the noise such as the power spectrum, 
the avalanche distributions and the pulse shape. In particular, the theory should be able to predict the values of 
the scaling exponents characterizing the distributions and possibly link the noise properties to microstructural details 
of the material under study. It is thus important to establish a controlled pathway from the microscopic reversal 
processes and the resulting macroscopic noise response. 

We can list here a set of crucial questions that need to be addressed by the theory: 

1. an explanation for the occurrence of scaling and quantitative estimates for the exponents. 

2. the origin of cutoff in the scaling regime. 

3. the effect of microstructural parameters (lattice structure, anisotropics, spin interactions) on the noise statistics. 

Clearly these questions are particularly broad and it is dubious that a single answer will be valid for all materials 
and experimental situations. The observation of scaling laws, however, can signify that relatively simple models 
could capture and reproduce a wide class of experiments. This is what is expected in critical phenomena, where 
universality implies that only the main symmetries and conservation laws influence the behavior of the system, while 
many other quantitative details are irrelevant. Thus, if the Barkhausen effect reflects the presence of an underlying 
critical phenomenon, one expects relatively broad universality classes spanning different materials, as it is indeed 
observed. 

Discussing the main theoretical approaches proposed to understand the Barkhausen effect, we will pay particular 
attention to the way each theory addresses the main questions reported above. In order to put the models into 
context, we first discuss the general properties of ferromagnetic materials in terms of the micromagnetic free energy. 
This represents a natural microscopic starting point to build up mesoscopic models and theories amenable to an 
analytic treatment and large scale numerical simulations. 

A. General properties of ferromagnetic systems 

A ferromagnetic materials can be represented as an ensemble of localized magnetic moments or spins, interacting 
between each other and with the external magnetic field H. The macroscopic properties of the material, such as 
the Barkhausen noise and the hysteresis loop, are due to microscopic rotations of the spins and could in principle 
be described by a microscopic theory. As a first step in this direction, we will discuss the interactions governing the 
dynamics of the local magnetization and the associated energetic contributions. In particular, we can write the energy 
of a ferromagnetic material as a sum of different terms 

E = Eex + E„i + Ean + Edis, (9) 

where E^x represents exchange interactions, E,n the magnetostatic energy, Ean the anisotropy and Edis the disorder. 
The detailed form of these terms, as a function of the local magnetization M{r) will be discussed in the following 
sections. 

1. Exchange energy 

The most important energetic contribution comes from exchange interactions, which are typically short-ranged and 
tend to align spins. For a set of spins s{fi) on the lattice, the exchange energy can be written as 

where J{x) decays rapidly for large x and the sum is over all atom pairs. Ea. llOl can be approximated in the continuum 
limit, through an elastic approximation which involves the removal of all small scale details. The set of spins s{ri) 
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are replaced by a continuous field M (r) and the exchange energy can be written as 

3 

E,,=A c^V^(VM„(f))^ (11) 
where A is the exchange coupling, which can be derived from J{x). 

2. Magnetostatic energy 

The magnetostatic energy is due to the interactions between the spins and the external field and to dipole-dipole 
interactions between different spins. For a uniformly magnetized sample the contribution of the external field to the 
magnetostatic energy is simply given by 

Em = ~^VM ■ H, (12) 

where M is the magnetization and V is the volume. 

In addition, we should consider the energy due to the demagnetizing field Hdm, the magnetostatic field generated 
by the local magnetization. In order to compute Hj^^j it is useful to define "magnetic charges" associated to the 
discontinuities of the normal component of the magnetization. For a surface separating two regions of magnetizations 
Ml and M2, the surface charge density is given by 

o- n • (All - M2) (13) 

where h is the vector normal to the surface. For instance, at the boundary of the sample, where the magnetization 
varies abruptly from Ms to zero, the charge density is given by Ms cos 9, where is the angle between the magnetization 
direction and fi. It is important to remark that magnetic charges are not physical charges, but are a convenient tool 
to compute the demagnetizing field and magnetostatic energy of magnetized bodies. Everything proceeds as in 
electrostatics provided we replace electric charge with magnetic charges and the electric field with the demagnetizing 
field which is thus given by 

i?.™(.-) = -V./^, (14) 

where the integral is over the surface separating two regions of constant magnetization. For a uniformly magnetized 
ellipsoid Hdm is constant and proportional to magnetization vector 

Hdm = ~kM, (15) 

where k is a geometry dependent demagnetizing factor. In this case the total magnetostatic energy is simply obtained 
by replacing H in Eq. ^1 hy H + Hdm ■ 

In general the demagnetizing field is not constant and Eq. E| has to be replaced by an integral 

^™ = -^J '^''^^^ • + ^*»(^))' (16) 

where Hdm is given by the contribution of surface charges (i.e. Eq. I14|l summed to the contribution due to the bulk 
variation of the magnetization. The bulk dipolar energy can be written explicitly as 

As in the case of the surface term Eq. 1171 can also be rewritten in terms of the density of magnetic volume charges 
p{f) = V • M(r) as 

^™ = ~8^./ • ^^^^ 
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3. Magnetocrystalline anisotropy and magnetoelastic energies 

The magnetization in a ferromagnetic material typically has preferential directions corresponding to the crystallo- 
graphic axis of the material. Therefore, it is easier to magnetize the sample along the easy directions. This observation 
can also be expressed in terms of the energy of magnetocrystalline anisotropy 




where Ma is the a component of the vector M and Kij is a symmetric tensor, describing the anisotropy of the 
material. In the simplest case of a uniaxial crystal, Eq.^| reduces to 

Ean^ J d^rKo{M ■ e)^ ^ j cPtKqM'^ sm^ (f) (20) 

where (j) is the angle between the easy axis e and the magnetization vector and Kq is the uniaxial anisotropy constant. 

The variations of the magnetization inside a ferromagnetic sample can cause deformation in the lattice structure, a 
phenomenon known as magnetostriction. Conversely, when an external mechanical stress is applied to the sample the 
magnetic structure can in principle be modified. To describe this effect it is useful to introduce the magnetoelastic 
energy, which in the most general form can be written as 




where aap is the stress tensor and Xap-yS is the magnetoelastic tensor. For a crystal with isotropic magnetostriction, 
under a uniaxial stress a the anisotropy energy takes the simple form of Eq. 1201 with Kq replaced by Kq + 3/2A(T, 
where A is the uniaxial magnetostriction constant. 

4- Disorder 

In the previous discussion we have considered only an homogenous system, in which the interactions are globally 
defined and do not depend on position. In general, however, different sources of inhomogeneities are found in virtually 
all ferromagnetic materials. The presence of structural disorder is essential to understand the fluctuations in the 
Barkhauscn noise, which would be strongly suppressed in a perfectly ordered system. The nature of the disorder can 
be inferred from the microscopic structure of the material under study. 

We can thus distinguish several contributions to the magnetic free energy due to the disorder: in a crystalline 
materials disorder is due to the presence of vacancies, dislocations or non-magnetic impurities. In polycrystalline 
materials we should add to these defects the presence of grain boundaries and variations of the anisotropy axis in 
different grains. Finally in amorphous alloys disorder is primarily due to internal stresses and the random arrangement 
of the atoms. It is important to notice that in the following we will consider the disorder as quenched (or "frozen"): 
it does not evolve on the timescales of the magnetization processes under study. It is not always simple to quantify 
the energetic contributions of the different sources of disorder but we can highlight here the main effects. 

The presence of randomly distributed non-magnetic inclusions give rise to a magnetostatic contribution, due to the 
magnetic charges that form at the boundaries of the inclusions [TQI . The energetic contribution due to a collection of 
those inclusions can be expressed as a local fluctuation of the dipolar coupling. A similar discussion can be repeated in 
the case of exchange interactions, leading to a fluctuating exchange coupling. The typical fluctuations of the coupling 
depend on the volume fraction v of the non magnetic inclusions: the typical strength of the (exchange or dipolar) 
coupling g{r) will be of the order of g ~ (1 — u)go and the fluctuations {{g — {g)Y) — vgg, where go is the coupling 
without impurities. This type of disorder is conventionally called of random-bond type and should also be present in 
amorphous alloys, due to the random arrangements of the atoms. 

In polycrystalline samples, each grain has a different crystalline anisotropy. In particular, the direction of the 
anisotropy axes will fluctuate in space and in the simple case of uniaxial anisotropics the anisotropy energy will be 
given by Ean = J (PrK{M ■ e(r))^, where e is a random fmiction of position reflecting the grain structure of the 
material under study. 

In magnetostrictive samples internal stresses play a similar role than anisotropics, as we discussed in the previous 
section. In particular, a random distribution of internal stresses produces an random energy of the type of Eg. 1201 with 
a random anisotropy constant K{r) cx Kq + 3/2AiT(r). While internal stresses can have multiple origins, it is possible 



23 



to compute explicitly their distribution in some particular cases, such as in a crystal with a random distribution 
of parallel dislocations. The stress in r due to a dislocation in the origin is given in cylindrical coordinates by 
(Jap = b^Cai3{0)/r, where b is the Burgers vector, fi is the shear modulus and Cap is an angular function depending 
on the dislocation type . The distribution of the internal stresses is formally given by 

Pi^cp) = / d^'^rDin, ...rN)Siaai3 - ^(Ta/3(r - n)), (22) 

k 

where D is the distribution of the position of the dislocations. The internal stress distribution has been evaluated 
in the case of a random short-range correlated dislocation arrangement: P{cr) is Gaussian for small stresses, with 
variance (Acr^) cx p, where p is the dislocation density and at larger stresses displays a power law tail, scaling as 
[Slf. It is interesting to remark that P{cr) is an experimentally accessible quantity, since it is directly related to the 
X-ray spectrum [s^ . 

We have discussed two common contributions to the disorder energy in ferromagnetic materials: random bonds 
and random anisotropics. To conclude the discussion, we remark that some diluted antifcrromagnets in a field can be 
effectively described by a random- field type disorder, with energy given by 

Edrs = Jd\hir)-M (23) 

where h is randomly distributed in space. Despite the fact that random field disorder is not present in ferromagnetic 
materials, it provides a useful theoretical ingredient to understand the role of disorder in hysteresis, as we will discuss 
in the following. 



5. Micromagnetic equations 

In the previous sections we have discussed the different energetic contributions to ferromagnetic materials. Collecting 
all the terms, considering for simplicity a uniaxial material, we can write the energy 

3 „ 

E=J2 d\[AiVMaf + K{Maeaf - {Ha + H'fl^)Ma], (24) 

a=l •' 

A, K and 6^ in principle depend on position and H^^^ is the component a of the demagnetizing field, discussed in 
section IrvTO 

The energy function reported in Eq. 1241 can in principle be used to compute the equilibrium properties of a fer- 
romagnetic materials. This problem is in general very complex and we do not attempt to treat it here. In fact, to 
describe the Barkhausen effect we are interested in the evolution of the magnetization in response of an increasing 
external field. The equation of motion for the magnetization can readily be obtained from Eq. 1241 

dM 

— =jMx H^ff (25) 

where 7 is the charge to mass ratio and H^ff = —5E/5M. Eg. 1251 does not include dissipation mechanism and would 
predict an indefinite precession of the magnetization vector. This problem can be overcome introducing phenomeno- 
logical laws for the dissipation. This results for instance in the Landau-Lifshitz-Gilbert equation whose numerical 
integration for different microstructures and boundary conditions is the subject of micromagnetism. 

It is tempting to use a micromagnetic approach describe collective effects arising in the Barkhausen effect from 
first principles. An attempt on this direction was made in Refs. 83, 84, 85] , where a Montecarlo version of the 
micromagnetic equation was used to simulate the Barkhausen effect in polycrystals. 

The model consists in a chain of N magnetic moments with orientations defined by two angles 0^ and Oi. The total 
energy of the system is given by 

E ~ ^ y^(sin^ ai + 2H cos 6i + msin^ 9i sin^ (pi — 2acos/?i,i+i), (26) 

i 

where is the angle between the magnetic moment and the local (random) easy axis, H is the magnetic field, a is 
the exchange-anisotropy energy ratio and m controls the magnetostatic energy. The last term takes into account the 
interaction between magnetic moments and /3ij is the angle between the moments in i and j. 
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The evolution of the magnetization is dictated by a Montecarlo algorithm. As the external field is slowly increased 
the magnetic moments rotate in bursts which are associated with the Barkhausen jumps. The results obtained for the 
jump distributions are in some range of parameters reminiscent of the experiments (i.e. a power law with a cutoff) but 
a quantitative comparison with experiments would be questionable. The model treats very accurately some aspects 
of ferromagnetic system, namely spin rotations, nucleation and random anisotropics, but completely neglects the 
fundamental aspects of the problem: dimensionality (the model is one dimensional), the range of interactions (only 
nearest- neighbor interactions are present, while in reality dipolar interactions give rise to long-range terms decaying 
as and demagnetizing effects. For these reasons, this model is not able to reproduce quantitatively the scaling 

properties of the Barkhausen effect. To this end, it would be necessary to study three dimensional models with 
long-range interactions, a task that promises to be computationally very demanding. 



B. Random energy models 

The difficulties arising in the direct analysis of the micromagnetic equations naturally lead to search for simplified 
models that could capture the essential mechanism underlying the Barkhausen noise without excessive complications. 
The main ingredient of any model for the Barkhausen effect is the presence of disorder, thus the simplest approach 
is to neglect most of the details of the micromagnetic energy, keeping only a uniform magnetostatic term Em and 
replacing the rest with a random function of the magnetization m 

E^F{m)-mH,ff, (27) 

where H^ff — H + Hdem- The rationale behind this approximation can be understood considering a single rigid 
domain wall dividing the sample in two domains. In this case the magnetization is proportional to the wall position x: 
m = Ms{2x/ L—1), where Ms is the saturation magnetization and L is the sample width. The sample is magnetized as 
the applied magnetic field pushes the domain wall across a disordered landscape represented by the random potential 
F{m). In order to define completely the problem, one needs to specify the statistical features of F{m), such as its 
distribution and correlations. 

Neel was the probably the first to use a random energy model to study the hysteresis loop properties at low fields 
in his theory of the Rayleigh loops [sgI Is^ . In the Neel model the random function was schematized as a series of 
parabolas with randomly distributed curvatures. The Barkhausen noise can be constructed from the random energy 
model by a visual inspection of the pinning field W{m) = — ^ as shown in Fig.|^ As the effective field is increased, 
the domain wall is pinned as long as i?e// < W{m), when Hf,ff reaches a local maximum of W(rn) the domain wall 
jumps forward until the condition Heff < W{m) is met. The Barkhausen sizes S are just the changes in magnetization 
Am occurring after these jumps or, in more mathematical terms, the first returns of the stochastic process described by 
VF(to). The random energy model in the form proposed by Neel, however, can not describe correctly the Barkhausen 
noise because the pinning field is essentially uncorrelated. Thus the first return distributions decays exponentially, 
while the experiments indicate a power law distribution. 

In order to obtain a power law distribution of returns the pinning field has to be long-range correlated. This point 
was first realized by Bertotti _8S,, ,-§9| , who included a Brownian pinning field in a random energy model. The model, 
widely known as ABBM |16Lri7| . yields a quite accurate description of the Barkhausen noise statistics. In the ABBM 
model the magnetization evolves according to an overdamped equation of motion 

dm , , , , 

— = ct- km + W{m), (28) 

where the external field increases at constant rate H — ct, k is the demagnetizing factor, and the damping coefficient 
has been set to unity rescaling the time units. The pinning field is a Brownian process where correlations grow as 

{{W{m)~W{m')f) = D\m-m'\. (29) 

The main predictions of the ABBM model can be obtained deriving Eq. 1)28(1 with respect to time and defining 
V = dm/ dt 

dv 

— = c — kv + vf{m), (30) 

where f(rn) = dW/dm is an uncorrelated random field with variance D. 
Expressing Eq. 1201 as a function of v and m only 

^^-~k + f{m), (31) 
dm V 
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we obtain a Langevin equation for a random walk in a confining potential, given by U{v) ^ kv — clog(w). Asymptot- 
ically, the statistics of v is ruled by the Boltzmann distribution 

P{v, m cxd) exp exp{~kv/ D). (32) 

The distribution in the time domain is obtained by a simple transformation and it is given by 

Piv) ^ P{v, t ^ oo) = ^c/.r(,%) ■ (33) 

A consequence of Eq. H33I) is that the domain wall average velocity is given by (w) — c/k. For cj D < 1 the velocity 
distribution in Eg. 1331 is a power law with an upper cutoff that diverges as fc — > 0. In this regime, the domain wall 
moves in avalanches whose size and durations are also distributed as power laws. For c/ D > 1 the motion is smoother 
with fluctuations that decrease as c/D increases. 

The avalanche size distribution equals the distribution of first return times of a random walk in the confining 
potential U{v). In the limit fc ^ 0, we can directly apply the exact result reported in Ref. ^9Q] for the first return 
times of a random walk in a logarithmic potential, which implies that the avalanche exponents depends on c and are 
given by 

T = 3/2-c/2L> a = 2 -c/D. (34) 

See also Ref. |H3 for alternative but not rigorous derivations of these results. 

The scaling of the cutoff of the avalanche distributions in the ABBM model can be obtained similarly in the limit 
c — !■ solving for the first return probability of a biased random walk and is given by 

Sq - (35) 
Using similar arguments one can also show that the cutoff of avalanche durations scales as Tq ~ [sof. 



C. Spin models 



In this section we discuss the attempts to understand the statistical properties of the Barkhausen effect by models of 
interacting spins, that can either be simulated on a computer or in some cases analyzed theoretically. This approach 
is conceptually simple and theoretically appealing, since it tries to derive directly macroscopic properties from a 
microscopic model, replacing the vectorial micromagnetic equations by simplified rules for the evolution of integer 
valued spins. The basis of this kind of approach is the notion of universality which is expected for critical phenomena: 
the values of the critical exponent and the shape of the scaling functions are independent on the microscopic details 
of the system, provided that the relevant symmetries of the problem are correctly taken into account. For this reason, 
we can expect to obtain the correct large scale behavior, such as the scaling properties of the Barkhausen effect, from 
relatively simple models. 



1. The random field Ising model 

Sethna et. al |9l| | have proposed the driven random-field Ising model as a prototype for hysteresis and avalanches 
at first-order phase transitions. The model was studied extensively analytically and numeric ally and the resulting 
scaling behavior was compared with the experimental results on the Barkhausen effect |HS ll^ lla, HE IS^I ■ 

In this models, a spin Si — ±1 is assigned to each site i of a d— dimensional lattice. The spins are coupled to their 
nearest- neighbors spins by a ferromagnetic interaction of strength J and to the external field H . In addition, to each 
site of the lattice it is associated a random field hi taken from a Gaussian probability distribution with variance R, 
P{h) = cxp{-h'^ /2R'^)/V2ttR. The Hamiltonian thus reads 

E^-J2Js^SJ-J2iH + h^)s^, (36) 

where the first sum is restricted to nearest-neighbors pairs. In the zero temperature dynamics used by Sethna et al. 
[9l| . the external field is ramped from — oo to cx) and the spins align with the local field (9^ 



Si = sign( J Sj + hi + H). 



(37) 
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FIG. 21: The Barkhausen jumps (denoted by arrows) from the perspective of the random energy modeL 

In this way a single spin flip can lead the neighboring spins to flip, eventually trigger an avalanche. 

For small values of the disorder R the first spin flips are likely to generate a big avalanche whose size is comparable 
to the system size, leading to a discontinuous magnetization reversal. On the other hand, a large disorder prevents 
the formation of large avalanches and the magnetization reversal is smooth. In fact, these two regimes are separated 
by a critical disorder Rc, where the avalanches are distributed as a power law. The behavior of the model close to Rc 
is very similar to that of equilibrium systems at the critical points and can thus be studied by standard methods, like 
mean-field theory and the renormalization group. 

A first qualitative picture of the behavior of the RFIM can be obtained by mean-field theory |7^|^. To this end, 
we consider Eq. 1361 extending the sums to all the pairs of sites and obtain 

E = ~^{JM + H + h,)si, (38) 

i 

where M = ^ ^ Si . The magnetization can then be obtained self-consistently as 

M = 1 - 2P{JM + H + h,<0) = l-2 / p{h)dh. (39) 

This equation has a single valued solution for R > Rc ~ ■\/(2/7r) J, while the solution becomes multivalued for R < Rc- 
In the second case, we have an hysteresis loop with a jump at He and a diverging slope dM/dH. For high disorder 
there is no hysteresis, but this is an artifact of mean-field theory that is not observed in low dimensions. At the critical 
point R = Rc the magnetization has a diverging slope, but no finite jump. In analogy with critical phenomena it is 
possible to obtain scaling laws describing the singularities occurring at the transition. In particular the magnetization 
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obeys the scaling law 

M^\r\^g±{h/\r\^'), (40) 

where r = {R — Rc)/Rc and h = {H — He) /He arc the reduced control parameters, and /3 = 1/2, S ~ 3. The scaling 
function g± is obtained as the root of a cubic equation. 

The distribution of avalanche sizes near the critical point can be computed analytically in the framework of mean- 
field theory. The result can be summarized in the scaling form 

P{S,r,h) = S-^V±{S/\r\-^/%h/\rn, (41) 

where the exponents take the mean- field values t — 3/2, a = 1/2. The scaling function 'P± can be evaluated exactly 
in terms of g± and is given by 

V±{x, y) = (1/V2^) exp(-x(l T Trg±{yf / /'^)- (42) 

In addition to the size distribution one can also define the distribution of durations T, which obeys an analogous 
scaling form 

D{T,r,h) ^T-'^V±[T/\r\-^'^'''''\h/\r\f^^), (43) 

where a = 2, is the correlation length exponent (^ ~ r^'^) and z is the dynamic exponent relating the correlation 
length to the characteristic time ~ Tq ). 

Mean-field theory has the advantage of being easily tractable and provides a qualitative picture of the behavior of the 
model, but the numerical values of the exponents are typically inaccurate in dimensions lower than the upper critical 
dimension, that for the RFIM it is equal to = 6. To overcome this problem, one can perform a renormalization 
group analysis with an expansion in e = 6 — d as in equilibrium critical phenomena. The renormalization group is 
quite involved and we do not discuss it here in detail, but just quote the main results (for a complete discussion the 
reader is referred to Ref. [t^)- The exponent determining the scaling of the order parameter have been computed 
to first order in e and are estimated to be (3 = 1/2 — e/6 and ^ = 3 -t- e. The avalanche exponent r displays only 
corrections to order and the cutoff exponent is estimated as ti = 1/2 — e/12. 

Large scale numerical simulations have been performed, to obtain reliable estimate for the critical exponents in 
three dimensions. The avalanche size distribution is most naturally computed integrating the avalanches over the 
entire hysteresis loop. The exponent Tint measured in this way can be related to the other exponents by a scaling 
relation Tint — t aj35. Typically the distribution displays an initial power law behavior and a cutoff that depends 
on the disorder R (see Fig. I22|l. The distribution for different values of R can then be collapsed using a scaling form 
of the type of Eq. ^] 69, 93- The best numerical estimate for the exponents (or exponent combinations) in three 
dimensions is r = 1.60 ± 0.06, r,„t = 2.03 ± 0.03, l/a = 4.2 ± 0.3, avz = 0.57 ± 0.03, a = 2.02, (3 = 0.035 ± 0.03, 

i/v = 0.71 m. 

Looking at the results above one may be tempted to perform a direct comparison with experiments. However, this 
should be done with a grain of salt. The RFIM model in fact predicts scaling only close to a critical disorder Rc 
while in general there is no reason to believe that this condition is fulfilled in experiments. To relate the Barkhausen 
scaling behavior observed in experiments to a disorder induced phase transition, such as the one discussed here, it 
would be necessary to show that the cutoff to the scaling behavior scales with the disorder strength. Unfortunately 
such a measure has not been performed yet. 

2. Random bonds, random anisotropics and other effects 

The RFIM represents probably the simplest model to study the effect of disorder in ferromagnetic hysteresis and 
Barkhausen noise. The model displays a disorder induced transition with associated critical exponents, which has 
been recently observed in experiments. While random fields can be a useful theoretical idealization, they are not 
present in ferromagnetic materials. Thus an important question to address is whether the basic phenomenology of 
the RFIM carries over to more realistic types of disorder. While the answer to this question appears to be positive, it 
is still controversial whether the disorder induced phase transition is universal with respect to the type of disorder. A 
second related question is whether the inclusion of more realistic ingredients in the RFIM results or not in a change 
of the scaling exponents. 

The first question was addressed by numerical simulations, considering random bonds, site dilution, and random 
(infinite) anisotropics. In the random bond Ising model (RBIM), the disorder enters into the exchange coupling Jij 
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FIG. 22: The integrated avalanche distribution in the RFIM in d = 3 simulated in a lattice of 320^ spins. The inset shows the 
data collapse. [From Fig. 1, pg. 4528] 



which is randomly distributed with mean Jq and standard deviation R. Hysteresis in the RBIM was first studied in 
the spin glass limit {Jq — 0) j9j|. In analogy with the RFIM, the RBIM displays a disorder induced phase transition 
in the hysteresis loop. Some of the scaling exponents, including those ruling the avalanche distributions, have been 
evaluated and their values differ slightly from those measured in the RFIM. For instance, some reported values are 
Tint — 2, z ~ 1.6 and ~ 1 |95j |. One should beware, however, that the lattice sizes used in this study are rather 
limited (up to 40^) compared to those used to estimate exponents in the RFIM (up to 320'^) _93J. 

In the site diluted Ising model (SDIM), the presence of non-magnetic impurities is treated by removing a fraction 
of spins in the lattice. The Hamiltonian is thus simply given by 

i? = -j ^ JijdCjSiSj - ^ Haai (44) 

ij i 

where = 0, 1 indicates weather the spin is magnetic or not. Typically, the sum in the exchange interaction is not 
restricted to the nearest neighbors but extends to a larger range. Simulations of the SDIM in d = 3 yield exponents 
compatible with those of the RFIM, namely Tint = 2.0 ± 0.2 |9q . 

The random (infinite) anisotropy model is obtained from a vectorial model (a lattice version of Eq. (|61|l 

E=-J2Js^■ s, - J2 K{s^ ■ + H ■ s„ (45) 

in the limit of ii' — > oo. Thus to minimize the anisotropy energy, in each site the spin should be parallel to the 
quenched random direction fii. The original vectorial model becomes then effectively of Ising type as can be seen by 
introducing a fictions spin ai — ±1 defining the sign of Si along its random axis (i.e. Si = fiiGi). In terms of the new 
spin variables the Hamiltonian can be rewritten as 

i? = - ^ Jijaiffj - ^ Htai, (46) 
{ij) ' 

with effective random bonds Jij = J hi ■ hj and random fields Hi = n^. The model above was simulated numerically 
with a zero temperature dynamics analogous to the one employed in the RFIM and with two different distributions 
for the anisotropy axis. The measured exponents are in good agreement with those measured in the RFIM. In 
particular for the avalanche distribution, the reported values are Tint = 2.06 ± 0.05 and Tint = 2.10 ± 0.05 depending 
on the distribution j93|. A recently proposed analysis, based on the renormalization group, of the vectorial random 
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anisotropy model shows that the model should be described generically by the RFIM fixed point [98j , thus supporting 
the conclusions of the simulations. 

From the discussion above, we conclude that the disorder induced critical point in the RFIM represents a broad 
universality class describing different sorts of spin models. This fact, if confirmed by larger scale simulations and 
renormalization group calculation, is particularly interesting because it allows to concentrate on a relatively simple 
model (the RFIM) to understand quantitatively a more generic scaling behavior. There is, however, an important 
point still to be considered: the effect of long-range dipolar interactions. As we discussed in section HV A 21 the effect 
is twofold, interactions between spins in the bulk and surface charges at the boundary of the sample. The bulk 
interaction has not been studied deeply in the framework of spin models, apart from some numerical simulations in 
two dimensions [99lllOClj| . From simple power counting considerations, we would expect that dipolar forces are relevant 
in the renormalization group sense, thus modifying the numerical values of the critical exponents. Therefore, unless 
one can neglect dipolar forces (this could be the case in hard rare earth alloys with strong anisotropics), the numerical 
results obtained in the RFIM would not apply to real materials. 

Surface effects have an even more drastic effect on the critical behavior discussed above. The presence of a demag- 
netizing field generates an effective correlation length, destroying scale invariance even at the critical point R — Re- 
Nevertheless, this effect is small, since the demagnetizing factor tends to vanish for very large system sizes. Sim- 
ulations of the RFIM in presence a demagnetizing field show remarkable features in the low disorder phase, where 
otherwise one expects the magnetization to change abruptly with a large jump. The demagnetizing field hinders the 
reversal process which splits in a series of avalanches of different sizes 101, 102]. The distribu tion of size is again a 
power law distribution and the cutoff is controlled by the demagnetizing factor (see Fig. |lOl| . It is important 
to observe that this power law scaling is not related with the disorder induced critical point, since this behavior is 
observed generically for R < i?c, but can be explained considering the dynamics of a domain wall. Indeed for R < Rc, 
at the coercive field a domain is nucleated and starts to grow. When the domain is large enough the demagnetizing 
field reduces the effective field and eventually the growth stops. At this moment the spins that will fiip when the 
external field is increased are likely to be those at the domain boundary. To understand this process in detailed one 
should thus analyze domain wall dynamics, as will be discussed in the next section. 

D. Domain wall models 

In order to minimize the exchange energy, a ferromagnetic material at sufficiently low temperature and zero field 
would display a uniform magnetization, oriented along an easy axis to lower the magnetocrystalline anisotropy en- 
ergy. In most cases, however, this arrangement is hindered by a very high cost in magnetostatic energy, due to 
the discontinuities of the magnetization at the sample edges. Thus it is energetically favorable to create domains, 
which in the simplest situation span the sample from end to end. In this condition, the magnetization process occurs 
prevalently b y do main wall motion, which can then be studied to recover the statistical properties of the Barkhausen 
noise |49l IHll IHI IT M lio l ITf)5| . 

1. Energetics of a domain wall 

We consider a fiexible 180° domain wall separating two regions with opposite magnetization directed along the 
z axis. In absence of disorder the domain wall would be fiat, but in general the wall will bend to accommodate 
pinning forces. If the surface has no overhangs, we can describe the position of the domain wall by a function h{f, t) 
of space and time (see Fig. I25|l . The energy associated with such a deformation field can be computed expanding 
the micromagnetic free energy discussed in Sec. IIV Al and we can thus split the energy into the sum of different 
contributions due to exchange and magnetocrystalline interactions, magnetostatic and dipolar fields, and disorder 

The contribution from the magnetostatic energy, if the external field H is applied along the z axis is given by 

Em - -2fioiH + Hdm)Ms J <fr h{r, t), (47) 

where we have included the demagnetizing field Hdm which is in general a complicated function of h(r, t), depending 
also from the sample shape. In the simplest approximation, however, the intensity of the demagnetizing field will be 
proportional to the total magnetization. Considering the field constant through the sample, it can thus be written as 

H,^^-^ j d\h{r,t) (48) 
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FIG. 23: The integrated avalanche distribution in the random infinite anisotropy model in d = 3 simulated in a lattice of 30'^ 
spins. The different curves correspond to different values of 6o, corresponding to the average angle between the field axis and 
the anisotropy axis [From Ref. Ha], Fig. 8, pg. 134431-6]. 



where the demagnetizing factor k takes into account the geometry of the domain structure and the shape of the 
sample and ^ is a volume factor. 

The interplay between magnetocrystallinc anisotropy and exchange interactions is responsible for the microscopic 
structure of the domain wall. While a very sharp change of the spin orientation has a high cost in exchange energy, 
a very smooth rotation of the spins between two domains is prevented by the magnetocrystalline anisotropy. The 
balance between these two contributions determines the width of the domain wall and its surface energy. The total 
energy due to these contributions is proportional to the area of the domain wall 



Ed^,, = 7- / d'rVl + |V/i(f,i)|2, (49) 



where 7^ ~ 2\/ AKa is the domain wall surface energy. Expanding this term for small gradients we obtain 

Ed^ = ^n,Sd^u + Y J d'r\Vh{r,t)\^, (50) 

where Sdw is the area of the undeformed wall. Thus the domain wall energy energy represents an elastic interaction 
that tends to keep the domain wall flat. 

An additional elastic interaction is due to dipolar forces, since the local distortions of the domain wall are associated 
to discontinuities in the normal component of the magnetization. As discussed in Sec. lIV A 21 this effect can be treated 
introducing magnetic charges on the domain wall surface (see Eq. I13|) . The surface charge density is zero when the 
magnetization is parallel to the wall and for small distortion it can be expressed as 

a(r) ^ 2Ms cos 9 ~ 2M, ^h^LlH (51) 

ox 



where 9 is the local angle between the vector normal to the surface and the magnetization. The energy associated 
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FIG. 24: Data collapse of avalanche size distribution P{S) ~ S f{S/So) simulated in the RFIM for R < Rc in presence of 
a demagnetizing field, using differe nt va lues of the demagnetizing factor k. The critical exponent r is given by 1.28, while 
So ~ k-'^^"" with l/dfc ~ 0.7 [From Fig. 1, pg. 233]. 



with this distribution of charges is given by 



27r|r — r 

which integrating by part can also be written as 



Ed= drdr — — (52) 



Ed = I d'rdPr'K{r- f'){h{f, t) - h{r' , t) f (53) 
where the non local-kernel is given by 



The interaction is long range and anisotropic, as can be seen by considering the Fourier transform 

K{p,q)^—^ / o (55) 

where p and q are the two components of the Fourier vector along z and y. In the preceding derivation we have 
implicitly assumed an infinitely strong anisotropy, so that the magnetization never deviates from the easy axis. In 
practice, however, the magnetization will rotate slightly from the easy axis, producing additional volume charges as 
discussed in section HV A 21 This effect leads to a slight modification of the interaction kernel in Eq. (|53|l 

K(p,q)^^ f (56) 

where Q = 1 + 2fioM^ / K is a material dependent constant |79l |. 

The disorder present in the material in the form of non-magnetic impurities, lattice dislocations or residual stresses is 
responsible for the deformation and the pinning of the domain wall. In general, disorder can be modelled introducing a 




FIG. 25: A domain wall separating two regions of opposite magnetization, 
magnetization across the domain wall produce magnetic charges. 



The discontinuities of the normal component of the 



random potential V^f, h), whose derivative gives the local pinning field r]{r, h) acting on the surface. In the particular 
case of point-like defects, such as non-magnetic impurities, the random force is given by 

V{r,h) = -^fp{f-fi,h- hi) (57) 

i 

where(ri, hi) are the coordinates of the pinning centers and fp{x) is the individual pinning force, which typically has 
a range comparable with the domain wall width Sw ~ y/ K/A. After coarse-graining at a scale larger than the typical 
distance between the pinning centers, this disorder can be replaced by a Gaussian random noise with correlations 

(ry(r, h)'q{r', h')) = 5^{r~ f')R{h - h') (58) 

where R{x) decays very rapidly for large values of the argument. It has been shown that the particular form of R{x) 
(i.e. due to random-bond or ra ndom-field types of disorder) has not a relevant effect on the scaling laws associated 
the domain wall dynamics 'l06]. 

Another possible source of pinning is due to local variations of the domain wall energy 7^, jj^ i due for instance to 
random anisotropics. Thus the domain wall energy becomes a function of position 7(r, h) which we can expand about 
its average + vi^j ^)- Introducing this expression in Eg. I5UI we obtain to lowest order an additional random term 
r]{ f, h), whose distribution and correlations can be directly related to the random anisotropics. 



2. Domain wall dynamics and depinning transition 



In most cases, the motion of the domain wall is strongly overdamped, since eddy currents cancel inertial effects. 
The equation of motion for the wall is thus given by 

^ dh{f,t) _ 6Ei{h{f,t)}) 

dt 6h{r,t) ' ^ ' 

where E{{h{'r,t)}) is the total energy functional, derived in the previous section, and T is an effective viscosity. We 
have neglected here thermal effects, since experiments suggest that these ar e no t relevant for the Barkhausen effect in 
bulk three dimensional samples [53. This could be different in thin films |l07j where thermal activated motion can 
be described by adding an extra noise term to the equation of motion. 
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Collecting all the energetic contributions, we obtain the equation of motion for the domain wall |50j. In order to 
avoid a cumbersome notation all the unnecessary factors can be absorbed in the definitions of the parameters. The 
equation then becomes 

^h^=H^kh + j^V'hir,t)+ (60) 



(fr'K{f- r'){h{r') - h{r)) + r]{r, h), (61) 

where the dipolar kernel K is reported Eq. (|54|l . the effective demagnetizing factor in given by fc = Afx^kM"^ /V and 
h = J d?r'h{f' ,t). Owing to the fact that the demagnetizing field term is just an approximation, the dependence of 
k on the sa mp l e shape and size can be quite complex . Var iants of Eq. H54(l have been extensively studied in the past 
|49ll5(ll5il7l[Tn ^ll0i[T0^ 

When the demagnetizing factor k is negligible, as for instance in a frame geometry, Eq. 1611 displays a depinning 
transition as a function of the applied field H and the domain wall moves only if the applied field overcomes a critical 
field He- For fields H > the domain wall moves with an average velocity v that scales as 

v^{H~ Hcfe{Hc - H), (62) 

where 9{x) is the step function. The critical behav ior associated to the depinning transition has been studied using 
renormalization group methods |l06l Ill4i lllfil Ill6j , which show that at large length scales the critical exponents take 
mean- field values ,49,. .,50.] . Th i s res ult is due to the linear dependence on the momentum of the interaction kernel 
(Eq. Isljl in Fourier space In general, if wc consider an interface whose interaction kernel in momentum 

space scales as K{q) = A^lql^, the upper critical dimension is given by = 2fi and the values of the exponents 
depend on fi. Notice that here we define d as the internal dimension of the domain wall. Hence in the present analysis, 
we have a two dimensional interface {d — 2) moving in a three dimensional medium. While in general we would expect 
mean-field behavior, there are situations in which the dipolar coupling could be neglected with respect to the domain 
wall tension 7^, where one effectively could observe the critical behavior associated with fi = 2. 

In general, for interfaces close to the depinning transition, the response to small variations of the applied field occurs 
by avalanches whose sizes S are distributed as 

P{S) ^ S-'f{S/So), (63) 
where the cutoff scales as Sq ~ {H — Hc)^^/" and is related to the correlation length ^ by 

^0 ~ (64) 

where C is the roughness exponent fFig. I26(l. The correlation length diverges at the depinning transition as 

^^{H-H,)-\ (65) 

which implies 

-^y[d+0- (66) 

The average avalanche size also diverges at the transition 

{S)^{H~H,)-\ (67) 

where 7 is related to r and a by 

a 

and to V and z by IlO(^ 

7=(1 + ^^C), (69) 

thus implying 

--2-4±4-. (70) 
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FIG. 26: The domain wall moves between two pinned configuration in an avalanche of size 



The other exponent relevant for the Barkhausen effect describes the distribution of avalanche durations 

F(r) ^ T-"5(T/To), (71) 

where the cutoff diverges at the transition as Tq ^ {H — Hc)^^/^. From Eq. (|65() and the relation Tq ^ we obtain 
A = l/zv and 



a — 1 



(72) 



Finally, an additional symmetry of the equation of motion can be used to reduce the number of independent exponents 
through the relation 



(73) 



The renormalization group analysis allows for a determination of the critical exponents in the framework of e = dc~d 
expansion. As discussed above, for fi — 1 we expect mean-field results r = 3/2 and a = 2, while for /i = 2 to first 
order in e one obtains t = 5/4 and a — 11/7. Recently a two lo op ex pansion has been carried out, allowing for an 
estimate to order O(e^) of the exponents (t = 1.24 and a — 1.51) |ll6j |. Finally, numerical simulations yield r = 1.27 
and a = 1.5 for fi = 2. 

The discussion above applies strictly to the case k = and scaling requires that H is close to He- If one would 
instead ramp the field slowly, sampling the avalanche distribution over all the values of the field, the result will be 
different. We need to integrate Eq. ESlover H {s^l 



Pint{S) 



dHS-^f{S{H - H^Y'") - S 



(74) 



1.72 for ^ — 2 and Tint = 2 for 
a + A. The numerical result is 



with Tint = T + (T. Using the values of r and a reported above, we obtain Tint - 
/i = 1. A similar discussion can be repeated for avalanche duration, yielding aint 
aint = 2.3 and for fi = 2 and aint = 3 for = 1. 

In a typical Barkhausen experiment, however, the situation changes since fc > and H is normally increased at 
constant rate. The effective field acting on the domain wall is given by H^ff — ct — kh. Thus, according to Eg. 1621 
when Heff < He the domain wall is pinned and the effective field increases. As soon as H^ff > He the domain wall 
acquires a finite velocity and h increases. As a consequence, the effective field will be reduced until the domain wall 
stops. This process keeps the domain wall close to depinning transition, but criticality is only reached in the limit 
c —> and A; ^ 0. A similar mechanism occurs in models of self-organized criticality |ll8| . 

The demagnetizing factor k represents a restoring force hindering the propagation of Barkhausen avalanches. Thus 
when /c > 0, we expect to find a finite avalanche size cutoff Sq. The domain wall can not propagate over distances 
larger than ^, defined as the length at which the interaction term is equal to the restoring force 



(75) 
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FIG. 27: The distribution of avalanche sizes in the short-range domain wall model at fc = when the field is swept from H = 
to H = H* < He- For H* He, we obtain Tint = t + a. 



which impUes 

^^{k/AKy" i^k = l/{f^ + d). (76) 

The avalanche size and duration distributions cutoff 5*0 and Tq can be obtained using the scaUng relations reported 
above 

So - D{k/Ak)-^/^\ 1/afe = (d + C)/(/i + d) (77) 

and similarly 

To^(fcMfe)-^^ Afc = z/(/i + d), (78) 



where D = ^ {rf) denotes the typ i cal fluctii ation of the disorder. Inserting in Eqs. H77I78(I the renormalization 
group results to first order in e |1061 

iniliT^ C = (2m - rf)/3 and z=_a - - 2d)/9, we obtain 1/ak = 2/3 and 
= (m~ (4/^^2(i)/9)/(/i-|-d). Using the results obtained to order jll6l |. one obtains instead for /i — 2, 1/crfe = 0.71 
and Afc = 0.345. The numerical values obtained from simulations arc in good agreement with this prediction (see 

Fig. 123 . 

The mean-field theory, which provides a good qualitative description of the depinning transition and describes 
quantitatively the data for d = dc (apart from logarithmic correctio ns), is obtained discretizing the equation of 
motion, coupling all the sites with the average domain wall position h '119\. The dynamics of such an infinite range 
model is described by 

^ = ct-kh + J{h-h,) + ij,{h), (79) 

where J is an effective coupling and i = 1, ...N. Summing over i both sides of Eq. H79|) . one obtains an equation for 
the total magnetization m 

d ^ 

^ = ct-km + J2v^{h)■ (80) 

i=l 
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FIG. 28: The cutoff 5*0 of the avalanche size obtained from numerical simulations for different demagnetizing factors. The two 
universality classes correspond to a kernel in the equation of the \q\'^ type in Fourier space, with jj, = 1 (long range dipolar 
interactions - up), and jj. — 2 (short range elastic interactions - down). 



This equation has the same form of the ABBM model (see Eq.^HJ 

To make the similarity more stringent, we should interpret ^^rji as an effective pinning W{m), with Brownian 
correlations. When the domain wall moves between two pinned configuration W changes as 

n 

W{m) - W{m) =^ Ar]„ (81) 
1=1 

where the sum is restricted to the n sites that have effectively moved (i.e. their disorder is changed). The total 
number of such sites scales as n '--^ Z'* and in mean- field theory is proportional to the avalanche size S = \m' — m\ 
(since S ~ Z"*"^^ and C = 0)- Assuming that the Arji are uncorrelated and have random signs, we obtain a Brownian 
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FIG. 29: (a) The distribution of avalanche sizes in the infinite-range model for different driving rates for A'^ — 32696. The fit 
of the power-law part yields r = 3/2 — c/2, with c = c/ D and 1/D ~ 1.2. (b) The corresponding distribution of avalanche 
durations. The power-law part is fit with an exponent a = 2 — c. 



effective pinning field 

{\W{m')-W{m)\^) = D\m' -m\, (82) 

where D quantifies the fluctuation in W. Thus the Brownian pinning flcld, introduced phenomenologically in the 
ABBM model, is an effective description of the disorder resulting from the collective motion of a flexible domain wall. 
As a consequence of this mapping, we expect that the infinite range model displays the sane frequency dependence of 
the exponents as in the ABBM model (Eq. I34|) . Simulations confirm this claim as shown in Fig. I^Hl 



V. DISCUSSION OF THEORETICAL RESULTS 



A. Avalanches distributions 



The models discussed in Sec. IIVI provide different result and theoretical interpretation for the pulse distributions. 
While we have reported already the results obtained for the various models, for the sake of clarity we summarize here 
again the main conclusions. We have distinguished between two general classes of models: spin models and domain 
walls model. Rigorously speaking the second class is contained in the first, since a domain structure is typically present 
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model 


r 


a 


avz 


1/a 


l/cTfe 




Tint 




RFIM d = 3 


1.60 ± 0.06 


2.05 ± 0.12 


0.57 ± 0.03 


4.2 ± 0.3 






2.03 


2.81 


MF (ABBM) 


3/2 


2 


1/2 


2 


2 


1 


2 


3 


LRDW 


1.5 


2 


0.5 


2 


0.65±0.05 


0.3±0.1 


2 


3 


SRDW 


1.27±0.02 


1.50 ±0.05 


0.57 ± 0.02 


2.2 ± 0.1 


0.72±0.03 


0.39±0.05 


1.72 


2.3 



TABLE II: The values of critical exponents estimated fr om s imulations of the RFIM [o^ and domain wall models with short- 
range (SRDW) and long-range (LRDW) interactions |2llll20(| . Mean-field (MF) results are the same for both classes of models 
and coincide with ABBM exponents. The exponent cr describes the scaling of the avalanche size cutoff with the reduced control 
parameter, R — Rc in the RFIM and H — He in domain wall models. For fe > and R < Rc the RFIM reduces to the SRDW 
model and those exponents apply. 

in spin systems, but the analysis of domain wall dynamics would be needlessly complicated. It is useful instead to 
treat separately the case in which domain wall propagation does not interfere with nuclcation. 

The two classes of models yield different interpretation of the origin of scaling in the Barkhausen pulse distributions. 
In spin models, scaling is ruled by the proximity to a disorder induced critical point. The cutoff to the power law 
scaling is determined by the distance from this critical point. The broadness of the critical region would justify the 
wide occurrence of scaling in experiments, without any apparent parameter tuning (69j . In domain wall models scaling 
is dictated by a different critical point: the depinning transition. The interplay between slow driving (c — > 0) and small 
demagnetization fc — > keeps the domain wall close to the depinning critical point. The cutoff of the distributions is 
then controlled by the value of fc, which in strip geometries is typically small. A collection of the numerical estimates 
of the Barkhausen noise distribution exponents for the RFIM and domain wall models is reported in Table UTI 

Depending on the problem at hand the appropriate model for the Barkhausen noise may vary. In soft magnetic 
ribbons the domain structure is typically controlled by dipolar interactions and the sample shape, imposing a parallel 
set of domain walls spanning the sample from end to end. In this condition, the depinning scenario should apply as it 
is testified by the experimental findings. The avalanche exponents are in good agreement with the prediction of the 
depinning model and the cutoff is controlled by the demagnetizing factor. 

In particular, experiments on polycrystalline SiFe are well described by the ABBM model or, equivalently, by 
a long-range depinning model. The exponents measured in amorphous alloys are instead well fitted by the short 
range interface model. This is due to the fact that in amorphous alloys the average anisotropy vanishes and thus 
dipolar interactions have no effect on the domain wall dynamics (i.e. Q in Eg. 1551 is very large). On the contrary, 
in polycrystalline materials, every crystal is the source of a relevant anisotropy and long range effects occur. This 
is confirmed by the fact that long-rEinge exponents are found when small crystallites are induced in an amorphous 
matrix by proper thermal annealing |2lj |. 

The scaling exponents measured in experiments on NiFe [sj] wires and Vitrovac 25] could be also explained in 
this framework if one assumes that the Barkhausen signal was not stationary. This could occur if the field is not 
swept up to saturation so that the effective field varies and it is not kept around He as in most cases. Then one 
would obtain integrated exponents Ti„t, ai„t instead of t and a. The striking similarity between the experimentally 
measured exponents and the theoretical ones (i.e. Tint = 1-72 and aint = 2.3) supports this conclusion. 

Spin models could instead apply to cases in which the domain structure is more intricate and the effect of dipolar 
forces negligible. A typical example would be an hard magnetic materials, with strong local anisotropies. Systematic 
investigation of the Barkhausen effect in hard magnetic materials has just been started and a quantitative comparison 
between experiments and theory is still not possible. 

An important problem to address, regardless of the theoretical framework, is the effect of the field rate on the 
Barkhausen noise statistics. As we discussed in section II V Bl the ABBM predicts a linear dependence from the 
exponents r and a from the field rate, encoded in the parameter c. This result applies as well for the mean-field 
domain wall model (see Fig. I29|l and agrees well with experiments on FeSi alloys. Simulations of the short-range 
domain wall model show that the scaling exponents do n ot ch ange continuously as in mean-field, but a peak appears 
for large values of the jump size and grows with the rate |l05j (see Fig. l3Uf) . While one could fit the result for higher 
driving rate with a different exponents, such a fit does not appear very reliable and a constant exponent coexisting 
with a growing peak seems a more appropriate approximation. A similar result is found in amorphous soft magnetic 
alloys, as shown in Fig. |S| 

The effect of the driving rate on the avalanche distributions was analyzed in more general terms in Ref. . A set 
of results can be derived under some basic assumptions: 

(i) The system is close to a critical point, so that in the adiabatic limit (c 0) the size and duration distributions are 
a power law. (ii) The avalanche signal is stationary under the time window considered, (iii) In the adiabatic limit, 
the average number of avalanche nucleation events per unit field increase is a smooth function of the driving field. 
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(iv) For low c the avalanche dynamics is independent of c. (v) The field increases between avalanche nucleations is 
independent of the avalanche sizes, (vi) The avalanche sizes are uncorrelated in time, (vii) Avalanches are nucleated 
randomly in space. 

All these assumption are reasonable both for the domain wall models and for spin models. Notice, however, that 
assumption (ii) restricts the derivation to a stationary signal and thus does not apply to distributions measured 
integrating along the entire hysteresis loop. 

Under the assumptions above, one can show that the effect of the driving rate depends on the value taken by the 
exponent a in the duration distribution for c —^ 0. For a > 2 the distributions are unaffected by the driving rate, 
while for a < 2 a peak appears at large sizes and durations. The interesting case is a = 2, which corresponds to 
a linear dependence of the exponents from the driving rate c, recovering the ABBM result in Eq. 1341 This general 
derivation provides an explanation of the results of experiments and numerical simulations, predicting the conditions 
under which drive dependent exponents should be observed. 



B. The power spectrum 



Explaining the power spectrum of the Barkhausen noise has been one of the main objectives of the past theoretical 
activity. In earlier approaches the powe r spe ctrum was considered as the result of a superposition of elementary 
independent events (see for instance [ml EH) without a clear relation to the microscopic magnetization processes. 
A step forward in this direction was made in the framework of the ABBM model . 

/^From Eq. 1281 one can derive a Fokker-Plank equation describing the evolution of the probability distribution 
P{v, t\vo) to have a velocity v at time t when the velocity at i = is wq E^ 

dP d A, , dvP\ , , 



Notice that, as in section llV Bl we have rescaled time setting a damping constant equal to one. The velocity correlation 
function is defined as 

G(i) = - c/A:)(«(0) - c/fc)) = / dvdv^iv ~ clk){vo - c/k)P{v, t\vo)P{vo), (84) 
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where P{vq) is given by Eq. Deriving Eg. 1841 and using Eq. we obtain an evolution equation for G{t) 

^--kG, (85) 

which can be easily solved yielding 

G{t) ^ c/kexp-kt. (86) 

The power spectrum has thus a Lorenzian shape 

This simple result contains several interesting predictions: (i) at large frequencies the spectrum decays as oj~'^. (ii) 
The scaling is cut off below a fr'equency loq = 1/k. (iii) The amplitude of the spectrum increases with the field driving 
rate c. These features describe qualitatively the shape of the experimentally measured power spectrum, but not 
quantitatively. For instance, the tail of the spectrum often decays with an exponent different than 2 and the cutoff 
frequency depends on k but not simply as 1/k. 

To reproduce the low frequency increase of the power spectrum, the ABBM model can be modified introducing a 
correlation length ^* in the pinning field 16]. This gives rise to a power spectrum of the type 

, ,2 



where Tc cx The ABBM predictions for the power spectrum can be compared with numerical simulations of the 
infinite-range domain wall model (Eq. I79|l for different values of c. The result in Fig. E] are, as expected in good 
agreement with the theory. 

Several attempts have been made in the past to link throu gh s cali ng relati ons the power spectrum exponent to the 
scaling exponents describing the avalanche distributions HE US 1^ IzE If03j | • The comparison between these results 



and experiments is, in most cases, not satisfactory. In any case, a simple exponent would just describe the tail of the 
power spectrum, while a complete theory should explain the entire shape observed experimentall y. T he pro blem was 
recently revisited in Ref. [74| where all the assumptions underlying previous scaling approaches |25l Isol l5lL l75l llOSj 
were carefully scrutinized and tested. Through a combination of analytical methods and numerical simulations, Kuntz 
and Sethna [2J| found a scaling relation describing the high frequency tail of the spectrum which seems to be in good 
agreement with experiments 65]. 

The derivation of Ref. (74] follows different steps involving some assumptions which can be tested numerically 
and experimentally. Close to criticality one expect that the average size for an avalanche of duration T scales as 
S{T) ~ yi/o'i^z and consistently the average shape scales as 

v{T,t) = T'/''''''-'Uape{t/T), (89) 

where v is the signal voltage, t is the time and fshape{t/T) is a (possibly universal) scaling function. 

The second assumption involves the probability P(v\S) that a voltage v occurs in an avalanche of size S. This 
probability scales as: 

P{v\S)=V-'f,oltage{vS^'''~^) (90) 

where fvoitage is another universal scaling function. Using this equation one can show that the energy E = {v^) of 
an avalanche of size s scales as E{S) ^ S^^'^"^. The analysis proceeds with the calculation of the voltage correlation 
function 

G(t) = / dt'{v{t + t')v{t)), (91) 



which can be decomposed into various contributions for each avalanche G(t\S). This quantity obeys a scaling law 

G{t\S) = S^'^^'-faitS-''''). (92) 
The contribution to the power spectrum is then obtained by a cosine transformation 



/"OC 

F{UJ\S)^ / dtCOs{cut)G{t\S) ^ S^fenergv{u;~^^^'''S). 

Jo 



(93) 
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FIG. 31: The power spectrum measured in the infinite-range domain wall model for different field rates c. The tail is well fitted 
by a Lorenzian function with amplitude proportional to c. The entire spectrum is in perfect agreement with the ABBM result. 



To obtain the power spectrum of the signal one should average Eg. 1931 over different avalanche sizes. A naive estimate 
of the integral resulting from the average of Eq. ESI using P{S) ^ S~'^ leads to the incorrect prediction p^L IHlL l7fil| 
F{uj) ^ tj-(3-'^)/<="'^^. The reason behind this incorrect result comes from the fact that fenergyix) ~ 1/a; for large 
X. Thus, if r < 2, the leading contribution to the averaged power spectrum comes from the upper cutoff in the size 
distribution, ultimately leading to the result 

F{iu) Lu'^^'"'\ (94) 

Using mean field exponents aiyz ~ 1/2, we recover the ABBM result F{uj) ^ lo^"^ . The result in Eg. 1941 agrees 
reasonably well with experiments and is confirmed by numerical simulations of the domain wall model (see Figs. 1311 
andE3- 



C. The avalanche shape 

The problem of the Barkhausen avalanche shape was first introduced in Ref. l74l and proposed as a strong test to 
characterize the universality classes of the Barkhausen noise. Subsequently, Ref. '60| analyzed more systematically 
the results of numerical simulations for the RFIM with and without demagnetizing field. The conclusion reached is 
that the average avalanche shape is in both case symmetric and close to a parabolic shape. The pulse shape changes 
slightly with and without demagnetizing field although a reasonable collapse of the two set of curves can be obtained 
multiplying one set by a factor (see Fig. I33() . Small differences emerge from a more detailed fitting procedure. A 
favorable comparison with experiments and a theoretical explanation of the numerical results is, however, still lacking. 

The gene ral problem of the shape of a sstochastic excursion was recently addressed in more generality by Baldassarri 
et al. |l23j | . The authors show that for a broad class of stochastic processes the average shape of the excursion is a 
semicircle. This occurs irrespective of the value of the scaling exponent (the equivalent of <jvz). Thus in this sense 
the shape would appear to be more universal than the exponents. On the other hand, the shape and symmetry of the 
average shape can be changed at will introducing correlations in the stochastic process |l23j and it is even possible 
to construct two stochastic processes with the sane scaling exponents, but different shapes. This calls into question 
the idea that the shape is a universal function |124| . 
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FIG. 32: The power spectrum measured in the short-range domain waU model for different field rates c. The tail is well fitted 
by a power law oj-^-\ in good agreement with Eg. 1941 which predicts an exponent 1/auz = 1.77. 




FIG. 33: The avalanche shapes for the RFIM with and without demagnetizing field (the second case is referred as front 
propagation). The shapes for different times have been collapsed using the value 1/zau = 1.72 and 1/zav — 1.75 in the two 
cases. An additional normalization factor allows to rescale the two set of curves together. [From 66], Fig. 1, pg. 046139-3] 
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The appr oach of Ref. |123| can be extended to compute the shape exactly of the Barkhausen pulses in the ABBM 
model |125) . Neglecting the contribution of the demagnetizing factor k, Eq. |^ describes the motion of a Id random 
walk in a logarithmic potential E{v) — —c log{v), where the magnetization s plays the role of time. The corresponding 
Fokker-Planck equation is 

aP(„,.)_c,^ c, A_ (95, 



ds dv \ V dv 

where P{v, s) is the probability to find the walk in v at s. We are interested in a solution of this equation with the 
initial condition P{v, 0) = S{v — vq) and an absorbing boundary at the origin v = 0. Following Ref. 90] we can express 
this solution in terms of modified Bessel functions. For the case < c < 1, which is the condition of the ABBM model 
to have power laws in the avalanche distribution (see Eq. ^ and Ref. 16]), the probability P{v, s\vq,0; c) for a walk 
starting at vq to be at v after a "time" s, in the limit uq — > is simply proportional to a power of v times a Gaussian 
with variance s, and the average excursion of the walk is then obtained as 



(v) cx ^s{S-.s) (96) 

This results means that the universal scaling function gshape is proportio nal t o a semicircle. 

It is also possible to calculate the function f shape in the time domain |l25j . By definition, the avalanche size s at 
time t is given by the integral of w(t, T) from time zero to time t: 

ft/T 

dt'v{t',T)(xT^ f{x)dx, (97) 



'0 Jo 
which provides an expression of s = s{t, T) as a function of t and T. If we now impose 

vit,T)^vis^s{t,T),S{T)) (98) 

we get an integral equation for f shape involving gshape 

fshape{x) (X gshape J g{x')dx'^ . (99) 

Using the form of gshape computed previously, we can solve this equation with the boundary conditions /(O) = /(I) — 
1. The solution is 

fshape{x) (X sin(7ra:) , (100) 

i.e. the universal function f shape is proportional to an arch of sinusoid. Finally for the normalized avalanche one 
obtains 



v{s,S)^S^-'''''tt^[s/S){1-s/S), (101) 
v{t, T) ^ T^/'"''~^t:/2 sin(7ri/T) . (102) 
D. Two dimensions 

The discussion presented so far focused on three dimensional systems, which is relevant for most Barkhausen noise 
experiments. The experimental results obtained in thin magnetic field require instead an analysis of two dimensional 
models. This task is not straightforward and several new features arise in two dimensions. The demagnetizing field, 
that plays a major role in determining the statistical properties of the Barkhausen noise, vanishes for very thin samples. 
For this reason, often we do not have have a simple parallel domain structure as in three dimensional ribbons and an 
analysis of single domain walls would not be appropriate. On the other hand, an approach in terms of spin models 
suffers for different difficulties, since there is no consensus on the presence of a disorder induced transition in two 
dimensions. 

The RFIM has been extensively simulated in d = 2, using extremely large systems sizes (up to (7000)^) [6^ 193^ . 
Despite this, it was not possible to obtain a reliable estimate of the threshold i?c for disorder induced transition. 
The simulations are consistent either with a vanishing threshold or with very large correlation critical exponents 
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(i.e. V = 5.3 ± 1.4, l/cr = 10 ± 2 with i?c = 0.54). An avalanche scahng exponent has also been estimated as 
Tint — 2.04 ± 0.04. Additional results and analysis have appeared in the literature, but are not particularly reliable, 
being based on much smaller system s izes. Typically the presence of a disorder induced transition is assumed and 
exponents are evaluated [60l Il26l Il27j . A similar discussion would apply other models analyzed in the literature, 
like the RBIM |9^, il2j|, the RFIM with site dilution 129] and dipolar interactions 99] . In summary, the large scale 
hysteretic behavior of disordered spin models in d = 2 is quite controversial and quantitative predictions seem at the 
moment not available. A possibility is that d = 2 is the lower critical dimensions, as it is expected for these models 
in equilibrium 130]. If this is the case, power law scaling may still be observed for low disorder, with the conjecture 
that r = 3/2 (6^. Nevertheless, we feel that at present the results of both experiments and models are still unreliable 
to a make a comparison. 

Even in thin films there are several cases in which a domain wall model may be appropriate. In particular, in the 
magnetoptical measurements of Ref. 's^j's^ one observes directly Barkhausen jumps due to domain wall motion. The 
equation of motion for a single domain wall in two dimensions is similar to the one in three dimensions (e. g. Eg. 161(1 
with some notable differences. For thin films with in plane magnetization the demagnetizing factor fc is virtually zero 
and in Fourier space the dipolar kernel (Eq. I54|l does not scale as q but as \o%(aq) where a is a small scale cutoff 
|l3ll fl32j . Apart from the logarithmic correction, the dipolar term is similar to the domain wall energy and thus 
in principle should not significantly affect the universality class of the depinning transition. Numerical simulations, 
however, yield slightly different results in the two cases (for instance C = 1 with long-range forces and C — 1.2 without 
|i33, 134]'). While the short -range interface model has been studied extensively in two dimensions, less is known for 
the long-range problem. A result that could be relevant for the Barkhausen effect would be the avalanche distribution 
exponent which has been numerically estimated for short-range models: r — 1.13 ± 0.02 *135^, r = 1.02 ± 0.29 ,136li] . 
The most reliable and precise estimate is probably r = 1.115 137]. 

Before comparing simulation data with the experiments one has to be certain of the quantity effectively measured. 
When the demagnetizing factor is zero, a ramp up of the field does not lead to a Barkhausen stationary signal with 
-ffe// — He, but the avalanches grow up to He- One would then measure an effective exponent Tint ^ t + a. Given 
the relation a = v / {d -\- C,) and C, = 1.25 and v — 1.33 115] we obtain Tint — 1-6. If instead of ramping up the field, 
one would instead keep it fixed just below He as in Rcf. 32, 33], then we would exp ect that the avalanche distribution 
scales with r. Notice, however, that the avalanche "size" S defined in Ref. |32.l33l] is not what we have defined, but is 
instead proportional to the change in the displacement A/i. Avalanches are only counted when the transverse length 
I (see Fig. I26f) is larger than the system size L. Hence S = LAft, with a constant L. Using scaling relation one can 
show that the exponent measured Th — ((C+ I)''' — 1)/C — 1-2, which is still lower than the estimated value but could 
be consistent with the data. 

Finally, we mention that spin models have been also used to study the dynamics of a single domain wal l inc luding 
overhangs. Simulations have shown a crossover from flat to fractal depending on the disorder strength |138| ]. The 
relation between these phenomenon and the Barkhausen noise still remains to be explored. 

VI. CONCLUSIONS AND PERSPECTIVES 

The Barkhausen effect has been studied for almost a century. While the essence of the phenomenon, i.e. an 
irregular magnetization reversal process occurring at the microscopic scale, was already indicated in the original paper 
by Barkhausen, a quantitative and detailed understanding was gained only recently. The key to this achievement 
is twofold: the injection in the field of methods and theories drawn from statistical mechanics of non-equilibrium 
critical phenomena and the identification of a standard experimental setup to collect reproducible and statistically 
significant data. These two aspects are tightly linked, since the necessity of rationalizing the experimental setup and 
to perform a detailed statistical analysis of the data was dictated by the theory, which was stimulated on its turn by 
the indication of scaling behavior coming from experiments. 

The successful understanding outlined above is restricted mostly to the Barkhausen noise measured in bulk soft 
magnetic alloys, ribbons in particular. These are particularly instructive from the theoretical point of view because of 
the relatively simple parallel domain structure, induced by the interplay between dipolar forces and the sample geom- 
etry. The magnetization process is mainly due to domain wall motion in a disordered landscape which is responsible 
for the avalanche dynamics leading to the observed noise. The scaling behavior describing the statistical properties of 
the noise (avalanche distributions, power spectrum) can be related to an underlying domain- wall depinning transition. 
Critical phenomena are typically associated to a certain degree of universality: the scaling exponents do not depend 
on the microscopic details of the system but only on general symmetries, dimensionality, and long-range properties 
of the interactions. One would thus expect to be able to classify materials in broad universality classes, sharing the 
same critical exponents, irrespective to the microscopic details, such as sample composition, grain size or disorder 
strength. 
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To some extent experimental results confirm this expectation: polycrystalline SiFe alloys with different grain sizes 
and Si content, together with partially crystallized Co based alloys were all found to share the same critical exponents. 
Amorphous alloys of different composition are instead grouped in a different class characterized by a different set of 
exponents. These two groups of exponents coincide with those associated to the depinning transition of a domain 
wall, with or without long-range dipolar interactions. The latter are enhanced by local anisotropics, induced by 
the polycrystalline grains and absent in the amorphous structure. The domain wall depinning scenario successfully 
predicts additional features of the Barkhausen effect, such as the dependence on the field rate and the demagnetizing 
factor, and can thus be used as a basis to understand the role of the microstructure on the macroscopic magnetization 
properties. 

The new frontier of the Barkhausen effect is represented mainly by thin films, as testified the recent literature. 
Inductive methods, that have played a leading role in noise measurements in bulk materials, become less effective 
as the thickness of the film is reduced to smaller and smaller sizes. Magneto optical techniques do not suffer of this 
problem and the improved resolution gained in time and space is allowing for precise noise measurements, joined to 
a direct visual inspection of domain wall dynamics. The effort in this direction, however, is still in process and a 
systematic analysis of the noise statistics in different materials has not yet been completed. Other techniques, such as 
magnetic force microscopy (139.| or Hall sensors 140] . are emerging as powerful tools to study the Barkhausen noise 
at the nanoscale. For instan ce, a recent experiment revealed Barkhausen jumps due to a domain wall depinning from 
the atomic Peierls potential |141| . As a note of warning, we note that these techniques reveal the magnetic properties 
of the surface and are thus not effective to study the crossover from bulk three-dimensional to thin two dimensional 
behavior. To this end, it would be desirable to improve inductive techniques, as recently proposed. 

Apart from the experimental problems outlined above, thin films pose new theoretical challenges. It is unlikely 
that a simple translation of three-dimensional models to two dimensions will provide a significant understanding of 
the Barkhausen effect in thin films. The domain structure is in most cases more complicated than in the bulk and 
a single domain wall model will be inappropriate. The general framework of non-equilibrium critical phenomena 
and the understanding gained in domain wall and spin models could nevertheless lead to a successful theoretical 
understanding. 

Finally, we would like to mention h ere hard ferrom agnetic materials for which few Barkhausen noise measurements 
have been reported in the literature |l42l Il43. Il44| . In these nucleation-type materials, the noise is not due to the 
motion of domain walls, but exclusively to grain reversal. Correlation between neighboring grains play a significant 
role 142], as well as the temperature. These systems, having large local anisotropics and nucleation fields, could be 
well represented by disordered spin models such as the RFIM. Recent studies of these models have contributed to 
renew the interest in the Barkhausen effect, but a clear experimental realization is still lacking. 
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APPENDIX A: TRANSLATION OF THE ORIGINAL PAPER BY BARKHAUSEN 

Two phenomena discovered with the help of the new amplifier 
by H. Barkhausen 

1. Noise during magnetization of Iron 

With the improvement of the new vacuum tube amplifier it is possible to generate easily a 10000 times larger 
current, which means a 100 million gain in power. With this apparatus, one can reveal electric or magnetic alternate 
fields that otherwise, due to their weakness, would be unaccessible to our knowledge. In some sense it is as if we 
discovered a microscope, and one with almost 10000 times magnification! Shortly after W. Shottky has indicated 
that, through a large amplification, one can, in some sense, hear the electrons flying in the amplifier tube, since the 
current they generate reveal the spontaneous fluctuations predicted by the kinetic gas theory. 

I discovered a similar phenomenon two years ago in occasion of an experiment that I did in collaboration with Dr. 
T. Tuczek. Iron produces a noise when magnetized. As the magnetomotive force is smoothly varied, the molecular 
magnets flip in jumps to their new position. Because of this, they generate irregular induction pulses in a coil wound 
around the sample, that can then be heard as a noise in a telephone. 
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FIG. 34: Figure 1 of the original paper 

I have now examined this phenomenon more closely. Fig. 1341 shows the experimental apparatus. The iron sample E 
to be studied was put inside a small coil S with a diameter of 25mm and 300 turns. One end of the coil was connected 
to a telephone T through a 10000 times amplifier V and the other end to a mirror galvanometer G. The latter was 
so strongly damped that a deflection was staying almost unaltered for several seconds. Then each modification of the 
amplitude in practice is proportional to the corresponding modification of the induction flux of the coil. Displacing 
the U-shaped 10cm long magnet M, one could observe at the same time the modification of the induction and the 
noise strength. I determined the contribution of the iron sample to the former, measuring and subtracting the smallest 
observable deflection of the galvanometer for a given displacement of the magnet, but in absence of the iron sample. 
The displacement was made by hand. For this reason, the precision was not so large. 

Next, it was shown that the noise became weaker for larger samples. 20mm large samples yielded effectively no 
more noise. As will be shown later, the reason for this is related to the too small variation of the induction (per unit 
area). The small magnet was not sufficient to magnetize strongly enough the large samples. So a pile of slabs gave 
rise to a weaker noise than a single slab coming from the same pile. The most suitable samples were found to be 1 to 
1/2 mm large wires or 5 to 10mm thin ribbons. Employing stronger magnets one could also use stronger iron samples. 

Later it was established on a large number of different samples that the noise was stronger as the iron was softer. 
Hardened steel gave almost no noise, while soft annealed iron yielded a particularly strong noise. The cause of 
this, however, was not the well known difficulty in magnetizing steel, that showed itself through a somewhat smaller 
deflection of the galvanometer given the same position of the magnet. The distinction in the noise existed also when 
soft iron was magnetized two to three times more weakly. Very soft iron still produced a noise when the magnet was 
moved at a l/2m distance. When it was placed in the close vicinity, the noise was so strong that it could be heard 
clearly in the telephone even without the amplifier. The distinction between different kinds of iron is so large that 
one could really build on this a method to investigate iron. 

The magnetic field variation was mostly done in such a way that the magnet was either displaced away by the side, 
turning it at 180° and then moving it back, or slowly turned away for 180° from a given position, in particularly close 
vicinity. The noise due to this motion was very strong in particular positions. The variations of the deflection of 
the galvanometer did not always correspond to the noise strength. When a variation in the direction R stopped at a 
point A (Fig. ESI, aU was suddenly quite and stayed perfectly quite as one made a variation back from A to B, even 
when the magnet was moved at will back and forth between A and B. As soon as the magnet overcome even slightly 
the position A, suddenly the full noise started again. In addition, by a further displacement in the direction R, the 
boundaries A and B were in some sense shifted away to the right. Then the boundary A remained again surprisingly 
sharp in correspondence to any motion between A and B. The boundary B was not sharp. The noise grew up slowly 
when one went back through B until C. At the same time, the boundary A became as well less sharp, in the sense 
sketched in Fig. 1361 The galvanometer also showed quite strong variations in the region between A and B. No jump 
could be detected during its motion, when A was crossed and the noise suddenly started. Clearly the observation 
method was not very precise. 

For thick samples it would be perhaps the same for a stronger field, but the region A-B was not reachable with the 
changes produced by the present magnet. Thus, as remarked above, with this no noise could be generated. 

One could maybe explain this peculiar behavior of iron in the following way: the individual molecular magnets group 
themselves into assemblies of various sizes with different degrees of stability. Noise results only from the disruption 
or formation of greater assemblies. For a motion corresponding to Fig. 1351 only small assemblies are modified in the 
region A-B and noise is does not start. For a motion across B in the direction -R, little by little also medium and 
large assemblies get involved. In the direction -t-R, all small and middle size assemblies have been already dissolved 
in the first change coming from far away up to A. Only the largest and more stable assemblies in this direction are 
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FIG. 35: Figure 2 of the original paper 
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FIG. 36: Figure 3 of the original paper 



still remaining. Those will be first dissolved by a further change beyond A. 

In good agreement with this idea, moving across B one hears a uniform, low sough, while a loud crackling noise 
starts across A. For very slow variations, for instance when the magnet is moved far away from the sample, one can 
clearly distinguish isolated snapping hits. It appears that in steel the molecular magnets do not have the ability to 
join in large assemblies. In this material one only hears the low sough. Also the behavior described in Fig. 1361 is in 
agreement with this idea. The weak assemblies, magnetized in the backward motion (from B to C in Fig. I35|l . are 
overcome again in the direction R, before A is reached. 

A phenomenon occurring in a stripe, thin as paper, of soft iron sheet is still non clearly understood. A turn of the 
magnet by 180° gave rise to the stronger noise when the magnet was placed at a distance of 10cm. By turning it at 
shorter distance, the noise was surprisingly very weak, even when the magnet was turned very slowly. 

A similar analysis would be desirable with the magnetodetector of Marconi and with the telegraph, also in view of 
practical applications. 

Dresden, Institut fiir Schwachstromtechnik, Mai 1919. 
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